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We study the massless field on D n = DP ^;Z 2 , where D C 
R 2 is a bounded domain with smooth boundary, with Hamiltonian 
H(h) — ~^2 x ^ y V(h(x) — h(y)). The interaction V is assumed to be 
symmetric and uniformly convex. This is a general model for a (2 + 
l)-dimensional effective interface where h represents the height. We 
take our boundary conditions to be a continuous perturbation of a 
macroscopic tilt: h(x) = nx ■ u + f(x) for x 6 dD n , u £ R 2 , and 
/ : R 2 — > R continuous. We prove that the fluctuations of linear 
functionals of h(x) about the tilt converge in the limit to a Gaussian 
free field on D, the standard Gaussian with respect to the weighted 
Dirichlet inner product (/, g)y = j D Yli PidifidiQi for some explicit 
/3 — f3(u). In a subsequent article, we will employ the tools developed 
here to resolve a conjecture of Sheffield that the zero contour lines of 
h are asymptotically described by SLE(4), a conformally invariant 
random curve. 

1. Introduction. The object of our study is the massless field on D n = 
D n ^Z 2 , with Hamiltonian U(h) = EfeeD* V(V/i(6)). Here, D C R 2 is a 
bounded domain with smooth boundary. The sum is over all edges in the 
induced subgraph of ^Z 2 with vertices in D and V/i(6) = h(y) — h(x) denotes 
the discrete gradient of h across the oriented bond b = (x,y). We take our 
boundary conditions to be a continuous perturbation of a macroscopic tilt: 
h(x) = nu ■ x + f(x) when x £ dD n for u G R 2 and /: R 2 — > R is a 
continuous function. We consider a general interaction V £ C 2 (R) which is 
assumed only to satisfy: 

1. V(x) = V(— x) (symmetry), 

2. < ay < V"(x) < Ay < oo (uniform convexity), and 

3. V" is L-Lipschitz. 

The purpose of the first condition is merely to simplify the notation since the 
symmetrization of a non-symmetric potential does not change its behavior. 
Note that we can assume without loss of generality that V(0) = 0. This 
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is the so-called Ginzburg- Landau effective interface ( GL ) model, also 
known as the anharmonic crystal. The variables h{x) represent the heights 
of a random surface which serves as a model for an interface separating two 
pure phases. 

The macroscopic behavior of the GL model has been the subject of much 
recent study. An important step in this development is the construction and 
classification of Gibbs states, which are infinite volume versions of the model. 
In two-dimensions, it turns out that the height variable h(x) diverges as the 
size of the domain tends to infinity, so in order to construct a Gibbs state 
one must first pass to the gradient field V/i. The existence and uniqueness of 
gradient Gibbs states is proved by Funaki and Spohn in [ 3], where they also 
study macroscopic dynamics. Deuschel, Giacomin, and Ioffe in [8] establish 
a large deviations principle for the surface shape with zero boundary condi- 
tions but in the presence of a chemical potential and Funaki and Sakagawa in 
[12] extend this result to the case of non-zero boundary conditions using the 
contraction principle. The behavior of the maximum is studied by Deuschel 
and Giacomin in [7] and by Deuschel and Nishikawa in [9] in the case of 
Langevin dynamics. The central limit theorem for linear functionals of the 
infinite gradient Gibbs states was proved first by Naddaf and Spencer for 
the static model with zero tilt in [21] and Giacomin, Olla, and Spohn handle 
the time-varying case with general tilt evolving under Langevin dynamics 
in [14]. 

The special situation in which the interaction is quadratic, i.e. V(x) = 
\x 2 , corresponds to the so-called discrete Gaussian free field (DGFF) or 
harmonic crystal. Many of the results in [7, 9, 12—14, 21] have been estab- 
lished separately for the quadratic case and often the results in this setting 
are more refined. The reason for the latter is that its Gaussian structure 
greatly simplifies its analysis; we will discuss this point in more detail later. 
For the DGFF, large deviations principles for the surface shape as well as 
a central limit theorem for the height variable was proved by Ben Arous 
and Deuschel in [1], the behavior of the maximum studied by Bolthausen, 
Deuschel, and Giacomin in [2] and by Daviaud in [4]. In a particularly im- 
pressive and difficult work, Schramm and Sheffield in [23] show that the 
macroscopic level sets are described by a family of conformally invariant 
random curves which are variants of SLE(4). 

Beyond having a Gaussian distribution, the main feature that makes the 
analysis of the DGFF tractable is that its mean and covariance are com- 
pletely described in terms of the harmonic measure and Green's function 
associated with a simple random walk. These objects are very well under- 
stood in the planar case [20], which often allows for very precise estimates. 
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The mean and covariance of the more general GL model also admit a repre- 
sentation in terms of a random walk (Helffer-Sjostrand representation or HS 
random walk [8, 15]). The general situation, however, is much more difficult 
because in addition to being non-Guassian, the corresponding representa- 
tions involve a random walk in a dynamic random environment whose be- 
havior depends non-trivially on the boundary data. The hypothesis that V 
is uniformly convex is helpful in assuring that the jump rates of the HS ran- 
dom walk are uniformly bounded from above and below. This implies that 
its Green's function is comparable to that of a simple random walk, which 
in turn allows for rough (up to multiplicative constants) variance estimates 
and, more generally, centered moments in terms of the corresponding mo- 
ments for the DGFF (Brascamp-Lieb inequalities). Furthermore, that the 
jump rates are bounded means that the Nash-Aronson and Nash continuity 
estimates apply, which give some rough control of the off-diagonal covariance 
structure. 

1.1. Main Results. The main result of this article is the following central 
limit theorem for linear functionals of the height. Let D C R 2 be a bounded 
domain with smooth boundary. For k > 0, let H K (D) be the Sobolev space 
of degree k and H~ K (D) its Banach space dual. 

Theorem 1.1 (Central Limit Theorem). Let f : R 2 — > R be a contin- 
uous function. Suppose that D n = D n ^Z 2 and that h n is distributed ac- 
cording to the GL model on D n with h n {x) = y> n {x) + fix) for all x E dD n 
where <p n {x) = nu ■ x. Let a u (b) = E[V"(r?(6))] where r] has the law of the 
Funaki-Spohn state with tilt u. Define the linear functional 

£v D (2) = E au(b)Vg(b)Vih n - <p n )(b) for g G H K (D). 

beD* n 

For every k > 4, the law of on H~ K {D) converges weakly with respect 
to the weak topology of H~ K (D) to a Gaussian free field on D, the standard 
Gaussian with respect to the weighted Dirichlet inner product (<7i,<?2)y = 
Id X^i Pidigxdig2, with boundary condition f, where iPi,^) is proportional 
to (E[V // (77(0,ei))],E[V"(r ? (0,e 2 ))]). Ln particular, p x = fi 2 > when u = 0. 

We will recall both the notion of a Sobolev space and the GFF in subsection 
2.2. We will also review the notion of a Funaki-Spohn state (as well as a new 
construction for zero-tilt) in subsection 4.2. 

As we mentioned earlier, central limit theorems have already been estab- 
lished by Naddaf and Spencer in [21] and Giacomin, Olla, and Spohn in [14] 
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for linear functionals of infinite gradient Gibbs states of the GL model. Both 
articles are based on the beautiful observation that the CLT can be reduced 
to a homogenization problem using the HS representation. The proof in [14] 
has more of a probabilistic flavor while the approach in [21] is to use PDE 
techniques. The reason that these results are restricted to infinite gradient 
Gibbs states is that the homogenization techniques they employ fail to carry 
over to the finite case. In particular, the main step in [14] is a proof that 
the macroscopic covariance structure of the gradient Gibbs state for the GL 
model is the same as that in the GFF by showing that the HS random 
walk converges in the limit to a Brownian motion. The key tool here is the 
so-called Kipnis-Varadhan method [19], which is to represent the random 
walk as an additive functional of the environment from the perspective of 
the walker. When the environment is an infinite, stationary, ergodic Markov 
process then it remains so when viewed from the walker, thus the conver- 
gence to Brownian motion is a consequence of Corollary 1.5 of [19]. If the 
environment is finite and, in particular, not ergodic with respect to shifts, 
this approach can no longer be used since the environment viewed from the 
particle is not ergodic. 

The covariance matrix of the limiting Gaussian in [ ] is given in terms 
of a complicated variational formula. It is therefore not explicit, except in 
the case of zero-tilt where it is possible to argue that it is proportional to 
the identity using rotational invariance. The proof of Theorem 1.1 gives 
the covariance matrix explicitly, up to a multiplicative constant, which is 
another new result for gradient Gibbs states. 

The careful reader may note that the convergence in Corollary 2.2. of 
[14] is in H~ K (D) for k > 3 while we require k > 4. The reason for the 
distinction is that we assume only continuity of the boundary condition /. 
This forces us to perform an extra integration by parts, which in turn puts 
an extra derivative on the test function. In the more restrictive setting of C 1 
boundary conditions, our proof also gives convergence in H~ K (D) for k > 3. 

The main step in our proof is motivated by the Markovian structure 
enjoyed by the quadratic case: the law of a DGFF on D n with boundary 
condition f n is equal in law to that of a zero-boundary DGFF on D n plus 
the discrete harmonic extension of f n to D n . This property is a higher di- 
mensional analog of the fact that a random walk Xt on Z conditioned to 
satisfy X tl = x\ and X t2 = %2 for t\ < t2 has the law of Y t + H t where 
Yt is a random walk on / = [ti , ^2] conditioned to vanish at t\,t% and H 
is the discrete harmonic function on / with boundary values H(t\) = 
and Ht 2 = Xt 2 . Our next theorem is a quantitative estimate of the degree 
to which this property approximately holds for the GL model. Although we 
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state it as our second theorem, it is the key step in the proof of Theorem 1.1 
and much of the article is dedicated to its proof. In order to give a precise 
statement of this result, we first need to setup some notation. 

Suppose that D C Z 2 is a bounded subset of diameter R > 0. Fix A > 
and let B^-(D) be the set of functions <f> : dD — > R satisfying max x6 go 1 4>(x)— 

u ■ x\ < AilogR) 1 . For r > 0, let D(r) = {x £ D : dist(x,3L>) > r}. With 
4> € B^-(D), let P^ denote the law of the GL model on D with boundary 

condition <j). In other words, P^, is the measure on functions h : D — > R with 
density 

J2 V(V(fcV *)(&))) 

b£D* / 

with respect to Lebesgue measure on R'^L Here, h \/ <f) is used to denote the 
function 




1.1) hV</>(x) 



h(x) if x € D, 
(b(x) if x G 3D. 



We will write 0-^(f{x)) to denote the set of functions g for which there 
exists a constant c-^ depending on A but independent of R so that \g{x)\ < 
c x |/(x)|.For/3 = (/3 1 ,/3 2 ), let 

(A^/)(x) = f 3 1 (f( x +e 1 )+f(x-e 1 )-2f(x))+f3 2 (f(x+e 2 )+f(x-e 2 )-2f(x)) 

where e\ = (1,0) and e 2 = (0, 1). Note that A^ is the usual discrete Lapla- 
cian for f3 = (1, 1). 

Theorem 1.2. Suppose that £ B^-(D) and (3 = f3(u) as in the 
statement of Theorem 1.1. There exists C,e,5 > depending only on V such 
that if r > CR 1 ~ e then the following holds. There exists a coupling (/i^, h^) 
o/P^jP^, such that if h: D(r) — > R solves the elliptic problem A^h = 
with /i|g£)( r ) = h = h^ — h^ then 

P[hyth inD{r)} =G\(ir 5 ). 
When u = 0, we can take (3 = (1, 1) so that A@ is the usual Laplacian. 

One of the main challenges in the analysis of the GL model is the lack 
of useful comparison inequalities for its mean. The difficulty is that the 
only explicit formula is given in terms of the annealed first exit distribution 
of the HS walk [8]. It is not possible to extract any sort of asymptotic 
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contiguity of this measure with respect to the harmonic measure of simple 
random walk using only that the HS walk jumps with bounded rates, which 
is all that is required to prove comparability of the corresponding Green's 
functions hence also of centered moments with DGFF. Indeed, examples 
have been worked out in the continuum setting of diffusions in which the 
two measures are absolutely singular and that the support of the former has 
a fractal structure. The situation is further complicated in the setting of the 
HS walk since in addition to being dynamic, its jump rates also depend on 
the boundary conditions, hence it seems difficult to rule out pathological 
behavior whenever the walk gets close to the boundary and the boundary 
conditions are rough. 

Applying Theorem 1.2 to the special case ip(x) = u ■ x gives the following 
estimate of the mean, which we believe to be sufficiently important that we 
state it as a separate theorem. 

Theorem 1.3. Suppose that tp G B^(D). There exists C, e,S > such 
that ifr> CR l ~ e and f3 = j3(u) as in the statement of Theorem 1.1 then 
the following holds. If h: D(r) — > R is the A@ -harmonic extension o/E^/i 
from dD{r) to D(r) then 

max \E^h(x) -h(x)\ = O^iT" 5 ). 

When u = 0, we can take (3 = (1, 1) so that h is harmonic with respect to 
the usual discrete Laplacian. 

It is worth pointing out that both of these theorems place no restrictions 
on the regularity of the boundary conditions ip, ip nor the regularity of 3D. 

1.2. Sequel. This article the first in a series of two and will be a pre- 
requisite for the second. In the sequel, we will make use of many of the 
estimates developed here in order to resolve a conjecture made by Sheffield 
(Problem 10.1.3 in [24]) that the macroscopic level lines of the GL model 
converge in the limit to SLE(4); the case of quadratic potentials is proved 
by Schramm and Sheffield in [ ■]. The two papers together are meant to be 
fairly self-contained. 

1.3. Outline. The remainder of the article is structured as follows. The 
second section is a short discussion of discrete and continuum Gaussian free 
fields. We chose to include the former part of this section since the special 
Markovian structure of the DGFF is the inspiration for Theorem 1.2 and also 
to serve as an illustration of the complications associated with non-quadratic 
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interaction. In the latter part, we provide a brief description of the GFF, 
the standard Gaussian law on Hq(D). A much more thorough introduction 
can be found in [2 >]. In Section 3, we will give a formal introduction to 
the GL model, its Langevin dynamics as well as the HS representation, 
and the Brascamp-Lieb inequalities. In Section 4, we will explain how the 
Langevin dynamics can be used to construct couplings of the GL model 
and prove an energy inequality for the discrete Dirichlet energy of such a 
coupling. This section is concluded with an equivalence of ensembles result: 
the Funaki-Spohn shift-ergodic gradient Gibbs state can be realized as an 
infinite volume limit of models on finite domains. In Section 6, we will prove 
Theorems 1.2 and 1.3 using an entropy estimate which is based on technical 
estimates from Section 5. Finally, Section 7 is relatively short and deduces 
the CLT from Theorem 1.2. We conclude the article with two appendices 
containing useful estimates on discrete harmonic functions and symmetric 
random walks. 

2. Gaussian Free Fields. In this section we will introduce the discrete 
and continuum Gaussian free fields (DGFF and GFF). The reason that 
we include a discussion of the latter separate from the general case of the 
GL model is to emphasize its special Markovian structure, which is the 
motivation behind the ideas used in Section 6. 

2.1. Discrete Gaussian Free Field. Suppose that G = (V U dV,E) is a 
finite, undirected, connected graph with distinguished subset dV ^ and 
edge weights ui > 0. The zero-boundary discrete Gaussian free field (DGFF) 
is the measure on functions h : V U dV — > R vanishing on dV with density 



with respect to Lebesgue measure. Here, h V has the same meaning as 
in (1.1) and Zq is a normalizing constant so that the above has unit mass. 
Equivalently, the DGFF is the standard Gaussian associated with the Hilbert 
space Hq(V) of real- valued functions h onV vanishing on dV with weighted 
Dirichlet inner product 



This means that the DGFF h can be thought of as a family of Gaussian 
random variables (h, /)y indexed by elements / £ Hq(V) with mean zero 




as)v = E w ( 6 ) v /( 6 ) v ^)- 
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and covariance 

(2.1) Cov((h,f)%,(h,g)%) = (f,g)%, f,g€H*(V). 

Although perhaps non-standard since our Hilbert space is finite dimensional, 
this representation is convenient since it allows for a simple derivation of the 
mean and covariance of h. Let A w : V — > R denote the discrete Laplacian 
on V, i.e. 

A-/(x) = J>(6)V/(6) 

bBx 

and let G u {x,y) = (A a ') -1 lr x \(y) be the discrete Green's function on V. 
Summation by parts gives that 

U,9% = -£ f(x)A"g(x) = Y1 A "f( x )9(x) for f,g€ H^V). 

x&V x&V 

Thus 

h(x) = (h,l {x} (-)) L 2 = -(h,G"(x,-)) v , 

hence 

Cov(h(x), h(y)) = (G"(x, •), G"(y, -))v = G w (x, y). 

Suppose that W C V. Then Hq(V) admits the orthogonal decomposition 
Hl(V) = Mi@Mb@Mo where Mi, Mb, Mo are the subspaces of H^(V) 
consisting of those functions that vanish on V \ W, are A^-harmonic off of 
dW, and vanish on W, respectively. It follows that we can write h = hi+hs+ 
ho with hi G Mi, hs G Mb, ho G Mo where hi, hs, ho are independent. 
This implies that the DGFF possesses the following Markov property: the 
law of h\w conditional on h\ v \ w is that of a zero boundary DGFF on W 
plus the A^-harmonic extension of h from dW to W. In particular, the 
conditional mean of h\w given fa|y\vF is the A^-harmonic extension of h\gw 
to W. 

More generally, if <p : dV — > R, the DGFF with boundary condition cf) is 
the measure on functions h: V — > R with h\gv = 4> with density 

^exp (-I^(&)(V(^V</>)(6)) 2 

That is, h has the law of a zero boundary DGFF on V plus the A^-harmonic 
extension of (j) from dV to V. 
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2.2. The Continuum Gaussian Free Field. The GFF is a 2-time dimen- 
sional analog of the Brownian motion. Just as the Brownian motion can 
be realized as the scaling limit of many random curve ensembles, the GFF 
arises as the scaling limit of a number of random surface ensembles [1, 14, 
18, 21, 22], as well as the model under consideration in this article. In this 
subsection, we will describe the basic properties of the GFF necessary for 
our analysis. Let D be a bounded domain in R 2 with smooth boundary 
and let C^(D) denote the set of C°° functions compactly supported in D. 
We begin with a short discussion of Sobolev spaces; the reader is referred 
to Chapter 5 of [10] or Chapter 4 of [27] for a more thorough introduc- 
tion. With No = {0,1,...} the non-negative integers, when / G Cq°(D) 
and a = {ai,a 2 ) G N 2 , we let D a f = d^d^f. For k G No we define the 
H k (D)-noi*m 

(2-2) ll/ll!r*(J»= E / \ Da fW\ 2 dx 

\a\<k JD 

where \a\ = a\ + a 2 - The Sobolev space Hq(D) is the Banach space closure 
of C$°(D) under || 

' \\ii k (D)- If s > is not necessarily an integer then Hq(D) 
can be constructed via the complex interpolation of Hq(D) = L 2 (D) and 
Hq(D) where k > s is any positive integer (see Chapter 4 section 2 of 
[27] for more on this construction and also Chapter 4 of [17] for more on 
interpolation). A consequence of this is that if T: C^°(D) — > C£j°(D) is a 
linear map continuous with respect to the L 2 (D) and H k (D) topologies then 
it is also continuous with respect to H S {D) for all < s < k. For s > 
we define H~ S {D) to be the Banach space dual of Hq(D) where the dual 
pairing of / G H~ S (D) and g G Hq(D) is given formally by the usual L 2 (D) 
inner product 

(f,g) = (f,g)L*(D) = f(x)g(x)dx. 

JD 

More generally, for any s G R the fP(Z))-topology can be constructed ex- 
plicitly via the inner product 

(2.3) (f,g) s = |(l-Ar/ 2 /.(l-A) s / 2 5 ; 

see the introduction of Chapter 4 of [ 27]. We are using the notation A for 
the Laplacian on R 2 to keep the notation consistent since elsewhere in the 
article A refers to the discrete Laplacian. Here, 



(1 - Aff = ^[(1 + e? + £ 2 2 )W)] for p G R 
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where 

■77(0 = / e-*-*f(x)dx 

is the Fourier transform of /. We will be most interested in the space Hq(D). 
Fix a positive definite 2x2 real matrix A. An application of the Poincare 
inequality (Chapter 4, Proposition 5.2) gives that the norm induced by the 
weighted Dirichlet inner product 

(f,9)v= [ VMf/Stffor f,geC?(D) 

JD ij 

is equivalent to || • ||_ffim)- This choice of inner product is particularly con- 
venient because it is invariant under precomposition by conformal transfor- 
mations when A is a multiple of the identity. 

The A-GFF h on D can be expressed formally as a random linear combi- 
nation of an (•, -, )y-orthonormal basis (/„) of Hq(D) 

h=} j a n f n 

n 

where (a n ) is an iid sequence of standard Gaussians. Although the sum 

defining h does not converge in Hq(D), for each e > it does converge almost 

surely in H~ € (D) ([25, Proposition 2.7] and the discussion thereafter). If 

, a 

/, <7 G C^°(D) then an integration by parts gives (/, g)y = — (/, A g). Here, 
A = VAV. Using this, we define 

(h, /) v = -(h, A A /) for / G C^{D). 

Observe that (h, /)y is a Gaussian random variable with mean zero and 
variance (/, /)y. Hence by polarization h induces a map C^°(D) —> Q, Q 
a Gaussian Hilbert space, that preserves the Dirichlet inner product. This 
map extends uniquely to Hq{D) which allows us to make sense of (h, f) v 
for all / £ Hq(D). We are careful to point out, however, that while (h, -)y 
is well-defined off of a set of measure zero as a linear functional on (D) 
this is not the case for general / € Hq(D). 

Suppose that W C D is a smooth, open set. Then there is a natural 
inclusion of H%(W) into Hq(D) given by the extension by value zero. If 
/ G C^{W) and g G Cq°(D) then as (f,g)§ = -(f,A A g) it is easy to 
see that Hq(D) admits the (•, -)y-orthogonal decomposition A4 ® Af where 

Ai = Hq(W) and M is the set of functions in Hq(D) that are A -harmonic 
on W. Thus we can write 



h = h w + h W c = 2^ a n f n + 2^ (5 n g n 

n n 
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where (a n ), (/3 n ) are independent hd sequences of standard Gaussians and 
(f n ), {9n) are ortho-normal bases of M. and A/", respectively. Observe that hw 
has the law of the GFF on W, h\yc the A -harmonic extension of to 
W, and h\y and h\yc are independent. We arrive at the following proposition: 

Proposition 2.1 (Markov Property). The conditional law ofh\w given 
^\d\w i s that of the A-GFF on W plus the A -harmonic extension of the 
restriction of h on dW to W . 

This proposition will be critical in the proof of Theorem 1.1. It also al- 
lows us to make sense of the A-GFF with non-zero boundary conditions: if 
/: dD — > R is a continuous function and F is its A -harmonic extension 
from 3D to D then the law of the ^4-GFF on D with boundary condition / 
is given by the law of F + h where h is a zero boundary ^4-GFF on D. 

3. The Ginzburg-Landau Model. The Ginzburg-Landau V</>-interface 
(GL) model is a general effective interface model first studied by Funaki and 
Spohn in [ L ] and Naddaf and Spencer in [21]. Suppose that G = (VUdV, E) 
is a finite, undirected, connected graph with a distinguished set of vertices 
dV. Let V G C 2 (R) satisfy: 

1. V(x) = V(— x) (symmetry), 

2. < ay < V"(x) < Ay < oo (uniform convexity), and 

3. V" is L-Lipschitz. 

The law of the GL model on V with potential function V and boundary 
condition if) : dV — > R is the measure on functions h: V — > R with h\gv = ip 
described by the density 



with respect to Lebesgue measure and h\/ tp is as in (1.1). 

3.1. Langevin Dynamics. Consider the stochastic differential system (SDS) 




(3.1) dhf(x) = J2V(V(hf\/ip)( y b))dt + V2dW t (x) for x € D 

tea; 
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Fig 1. A typical realization of the GL model with zero boundary conditions on {0, . . . , 30} 2 
and potential function V (x) —x 2 +cos(x), sampled using a discretization of the SDS (3.1) 

where Wt(x), x S V, is a family of independent standard Brownian motions. 
The generator for (3.1) is given by 

where 

b&E 

is the Hamiltonian for the GL model. Thus it is easy to see that C is self- 
adjoint in the space L 2 (e - ^^), hence the dynamics (3.1) are reversible with 
respect to the law of the GL model. These are the Langevin dynamics. 

3.2. The Helffer-Sjostrand Representation. We showed in subsection 2.1 
that if V(x) = ^x 2 then the mean height is harmonic and that the covariance 
of heights is described by the discrete Green's function. Both of these quan- 
tities admit simple probabilistic representations: if Xt is a continuous-time 
random walk (CTRW) on G that jumps with uniform rate 1 equally to its 
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neighbors and r is the time it first hits dV, then 

Eh{x) = E x h(X T ) and Cov(h(x),h(y)) = E X [ l {Xs=y }ds 

Jo 

where the subscript x indicates Xq = x. The idea of the Helffer-Sjostrand 
(HS) representation, originally developed in [15] and reworked probabilisti- 
cally in [8, 14], is to give an expression for the corresponding quantities for 
the GL model in terms of the first exit distribution and occupation time of 
another CTRW. In contrast to the the quadratic case, the CTRW is rather 
complicated for non-quadratic V as its jump rates are not only random, 
but additionally are time varying and depend on the boundary data. Nev- 
ertheless, the HS representation is a rather useful analytical tool due to 
comparison inequalities (Brascamp-Lieb and Nash- Aronson) . 

Specifically, let hf solve (3.1) with boundary condition ifi. Conditional on 
the realization of the trajectory of the time- varying gradient field (V/if (6) : 
b £ D*), we let xf be the Markov process on G with time- varying jump 
rates V"(Vhf(b)). Let r = inf{t > : xf £ dV}. Let P, denote the joint 
law of (hf,xf) given Xq = x and E x the expectation under P x . 

Lemma 3.1. The mean and covariances of h^ admit the representation 

(3.2) Cov(^(x), h%)) = E x £ l {xf=y} ds 

(3.3) Eh^(x) = f 1 E x i;(X r T ^)dr. 

Jo 

We refer the reader to [8, Section 2] for a proof and also a much more 
detailed discussion on the HS representation. 

3.3. Brascamp-Lieb Inequalities. For v 6 RJ^I, we let 

<^> = J>(*)/^(x). 

The following inequalities, first proved in [3] and redeveloped probabilisti- 
cally in [8] , bound from above the centered moments of h^ with those of h* , 
where h* is a zero-boundary DGFF on G. Recall that ay,^4v are positive, 
finite constants so that ay < V" < Ay. 
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Lemma 3.2 (Brascamp-Lieb inequalities). There exists a constant C > 
depending only on ay, Ay such that the following inequalities hold: 

(3.4) Var((z/, h^)) < CYax({v, h*)), 

(3.5) Eexp((z/,^) -E(i/,/^)) < Eexp(G>, 

/or a// i/ G R |y| . 

We again refer the reader to [8, Section 2] for a proof. The Brascamp-Lieb 
inequalities allow for the following bound on the moments of the maximum 
which we will make use of many times throughout the rest of the article. 

Lemma 3.3 (Moments of the Maximum). Suppose that F C Z 2 is bounded 
and connected with R = diam(F). Let ( G B^-(F), /i** ~ P|,, and M = 
max xe i? \h^(x) — u ■ x\. For every e > and p > 1 we have that 

(EM*>) 1/p = Xp (R*). 

Proof. We may assume without loss of generality that A > 1. Combining 
the Brascamp-Lieb and Chebyshev inequalities, we have the tail bound 

P[M >t] = G\(exp(G\((logi?) X ) - t)). 

Furthermore, we have 

BM p < E|/i c (x)| p < O x (R 2 (\ogR) pX ). 

Consequently, 

EAf p < R pe + EM p l {M > ff} < R pe + (EM 2p ) 1 / 2 (P[M > R e }) 1 ' 2 
< R pe + O x ( J R(logi?) pS exp(c p - iiT)). 

Therefore (EMP)Vp = 0^ p {R € ), as desired. □ 

4. Dynamics. We now specialize to the case where G is a bounded, 
connected subgraph of Z 2 . We will write D for its vertices, dD = {x G Z 2 : 
dist(x,D) = 1} for its boundary, and D* = {b = (xb,yb) G (Z 2 )* : G 
-D} for its edges, where (Z 2 )* denotes the set of edges of Z 2 . Finally, let dD* 
be the set of edges that are either contained in dD or intersect both dD and 
D. The Langevin dynamics are extremely useful for constructing couplings 
of instances of the GL model with either different boundary conditions, 
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defined on different (though overlapping) domains, or both. Suppose that 
h^,h^ are solutions of (3.1) driven by the same Brownian motions with 
boundary conditions tp,ip, respectively. Let ip = ip — ifi and h = — . 
Observe that 

(4.1) dht(x) = 5>'(V(fc? V ^){b)) - V'(V(hf V $)(b))]dt 

b3x 

Let 

(4.2) c t (b)= f V"(V(hf + sh t )(b))ds and £ t f(x) = Y / c t (b)Vf(b). 

Then we can rewrite (4.1) more concisely as 

(4.3) dh t (x) = £ t h t (x)dt. 
The following is [13, Lemma 2.3]: 

Lemma 4.1 (Energy Inequality). For every T > we have 



x&D J ° b&D* 

(4.4) <C[^|SoWP+ [ T £ ®(xbWht(b)\dt) 

\xeD ^° bedD* J 

for C > depending only on ay , Ay . 

More generally, if ft solves dtf = Ctft then f t also satisfies (4.4). 

4.1. Coupling Bounds. The purpose of the next lemma is to show that 
limT_ ! . 00 (/j^, hip) gives the unique invariant measure of the Markov process 

(hf , hf), i.e. where hf, hf both solve (3.1) with the same driving Brownian 
motions. 



Lemma 4.2. 

1. The SDS (3.1) is ergodic. 

2. More generally, any finite collection h 1 , . . . , h n satisfying the SDS (3.1) 
and driven by the same family of Brownian motions is ergodic. 
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Proof. Part (1) follows immediately from Lemma 4.1. Indeed, the Poincare 
inequality implies that there exists cd > such that for all functions 
/ : D — > R with f\so = we have 

x£D b£D* 

Thus Lemma 4.1 implies that if (hf,hf) both solve (3.1) with the same 
Brownian motions and boundary data (though with possibility different ini- 
tial distributions) then with ht = hf — hf we have 

£(M*)) 2 + j- f Y.^(x)fdt <cY^ (Mb)) 2 . 

ieD D Jo x&D xeD* 

In particular, the integral is bounded as T — > oo which implies 



/•oo 

lim / V (h t (x)) 2 dt = 0. 



" 1 xeD 
The energy inequality also gives 

i^(M*)) 2 <E^^)) 2 

for all < t < T. Therefore ht as t — > oo. This proves (1). 

We will now prove part (2). In the interest of keeping the notation simple, 
we will prove the result in the special case n = 2. Suppose that (h\,h 2 ) 
solve the SDS (3.1) driven by the same Brownian motions with boundary 
conditions ijj l ,ijj 2 but with arbitrary initial conditions (/iq,/iq). We know 

that h\ converges in distribution P^, by part (1). Consequently, the pair 
(h l t ,h 2 ) is tight. 

We will now prove the existence of the limit linit_ 5 . 00 (/4 , h 2 ). Suppose that 
(Tfc) and (S^) are arbitrary increasing sequences diverging to infinity. Fix k 
and assume that < Sk- Let (h],h 2 ) solve (3.1) where h l = h' l Sk __ Tk . Let 

(hj,h 2 ) be another pair of solutions to (3.1) with the same boundary and 
initial conditions of {h},h 2 ) but driven by the same Brownian motions as 

(h l ,h 2 ). Then (/%>^T fc ) = {h\ k ,h 2 Sk ) and C h r k ^T k ) = ( h T h i h T h )- % the 
energy inequality, 

1 [ Tk £ \vChl-K)(b)\ 2 dt<^Y,\^( x )-K( x )\ 2 

JO ,^ r>* ± k ^ 



<^E(i^wi 2 + i^-T fc (x)i 2 ). 
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As Tfc — > oo, the first term in the summation on the right hand side clearly 
converges to zero almost surely. The second term in the summation converges 
to zero in distribution since h\ is tight. Consequently, for every 5 > there 
exists /co sufficiently large so that for all k > ko with Sk > Tk we have 



1 



.. Y, \V(hl-ht)(b)\ 2 dt>5 
*k Jo b&D * 



< 5. 



Let e > be arbitrary and fix 5 = e/co so that if / : D — > R is an arbitrary 
function vanishing on dD with X^beD* \^f(b)\ 2 — <5 then YlxeD \f( x )\ 2 — e - 
Assume that k is sufficiently large so that with probability 1 — e we have 



1 



b&D* 



\V(hl-hi)(b)\ 2 dt < 5. 



Then there exists (random) t G [0,T k /2] such that J2b&D* \^(Ko~ h t )( b )\ 2 ^ 
5 hence X^eD l^loO^) ~~ K ( x )\ 2 — e with probability 1 — e. Applying the 
energy inequality once again yields 

£ ^(x) -tf Tk {x)\ 2 <C^\hl (x)-hl (x)\ 2 < Ce 

xeD x£D 

with probability 1 — e. Of course, we can do the same if > Sk- Therefore 

we conclude that the subsequential limits of (hj, h 2 ) are unique, hence \i = 
lim^oo^, h 2 ) exists. The same argument also implies that for any s > we 
have limt_>oo(^<s+t> ^s+t) exists and has the same distribution as fi. Therefore 
H is a stationary measure. 

To finish proving the lemma, we just need to establish uniqueness. Sup- 
pose that each of the pairs {h\,h 2 ),{h\,h 2 ) solve the SDS (3.1), h\,h\ all 
driven by the same Brownian motions. Suppose further that both pairs have 
stationary initial conditions. Then we can use the energy inequality exactly 
in the same manner as in the proof of part (1) to deduce that \h\ — h\\ -4-0 
as t — > oo. Since (hj,h 2 ), (hj,h 2 ) are stationary, it therefore follows that 

{hlh 2 ) = Chlhl). □ 

Suppose that (hf , hf) is the stationary coupling of instances of the model 
with boundary conditions i/j,tp 6 B^-(D), respectively. Letting h = — h^, 
the Caccioppoli inequality (B.4) implies that 

/2r 2 si r>2r' 2 

J2 \Vh t {b)\ 2 dt<~ J2 M*)\ 2 dt 

2 b£B*(x ,r) x£B(x ,r) 
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for r > and xq 6 D with B(xq, 2r) C D. The maximum principle implies 
that h attains its maximum on dD, hence h = (3^((log i?) A ). Consequently, 
taking expectations of both sides of (4.5) and using the stationarity of the 
dynamics yields 

£ E|V^(6)| 2 <O x ((lo gj R) 2 *). 

beB*{x ,r) 

Recall that D(r) = {x £ D : dist(x,9-D) > r} and let e > be arbitrary. 
Then D(R 1 - e ) can be covered by 0(R 2 /R 2 - 2e ) = 0(R 2t ) balls of radius 
R}~ e , all of which are contained in D. Therefore 

£ E|v^(6)| 2 < o^ri*). 

beD*(R 1 ~ i ) 

Let D(R 1 ,R 2 ) = {x £ D : Ri < dist(x,dD) < R 2 }, r k = kR 1 '^, and 
Dfc = D{R 1 ~ <L + r-fc, + rfc + i). Note that we can write 

R ic -l 

X) E|V^(6)| 2 = E E E |V^(b)| 2 = 0a(^)- 
This implies there exists < k < i? 4e — 1 such that 

£ e|vm&)| 2 < Oy(ir e ) 

We have proven: 

Lemma 4.3. Suppose that (h^,h^) is a stationary coupling of two solu- 
tions of the SDS (3.1) driven by the same Brownian motions with ip,^ G 
B^-(.D). For every e > f/iere exists i? 1_e < i?i < 2R 1 ~ € such that with 
R2 = R\ + -R 1 ^ 56 we have that 

B\Vh t (b)\ 2 = 0^(R- e ) and ^ E|V/^(o)| 2 =0^(R 3e ). 

beD*(R 1 ,R 2 ) 6eD*(i?i) 

This lemma will be particularly useful for us in Section 6 in order to con- 
struct an intermediate coupling of , P^, exhibiting pointwise regularity 
near dD with high probability. 
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4.2. Gradient Gibbs States and Equivalence of Ensembles. By the re- 
verse Brascamp-Lieb inequality [8, Lemma 2.8], it follows that if D n is any 
sequence of domains tending locally to the infinite lattice Z 2 and, for each n, 
h n is an instance of the GL model on D n then Var(/i n (x)) — > oo as n — > oo. 
This holds regardless of the choice of boundary conditions, which suggests 
that it is not possible to take an infinite volume limit of the height field h n (x). 
However, the Brascamp-Lieb inequality (Lemma 3.2) gives that Var(V/i n (6)) 
remains uniformly bounded as n — > oo, indicating that it should be possible 
to take an infinite volume limit of the gradient field. 

Working with gradient rather than height fields, though unnecessary for 
d > 3, is convenient since it allows for a unified treatment of Gibbs states 
for all dimensions. Let X be the set of functions n: (Z d )* — > R. Let J- = 
o~(i](b) : b G (Z d )*) be the a-algebra on X generated by the evaluation maps 
and, for each D* C (Z d )*, let F D * = a(r](b) : b G D*) be the a-algebra 
generated by the evaluation maps in D*. Suppose that D C Z rf is bounded, 
ip G X, and V(j) = <p. If h is distributed according to P^, then the gradient 
field V/i induces a measure P^,, on functions D* — > R. We call P^» the law 
of the GL model on D* with Neumann boundary conditions (p. Let [x be a 
measure on X and suppose that r\ has the law fi. We say that n is a gradient 
Gibbs state associated with the potential V if for every finite D* C (Z 2 )*, 



Fix a vector u G R^. A gradient Gibbs state /i is said to have tilt u if 
E M ry(x + bi) = u ■ e% where b% = (0, e^), e% a generator of Z rf , and x G Z d is 
arbitrary. We say that \x is shift invariant if \x o t~ 1 = /x for every x G Z d 
where t x : Zi d — > Z d is translation by x. Finally, a shift-invariant [i is said 
to be shift-ergodic if whenever / is a shift-invariant ^-measurable function, 
then / is //-almost surely constant. 

Funaki and Spohn in [13] proved that the shift-ergodic Gibbs states are 
parameterized according to their tilt u; from now on we will refer to the law of 
such as SEGGS U . The natural construction is to take an infinite volume limit of 
gradient measures P^* as D* tends locally to Z d with boundary conditions 
ip(b) = u ■ (yb — Xb) [ , Remark 4.3]. The difficulties with this approach 
are that P^,* is itself not shift-invariant and it is not clear that the mean 
gradient field approximately has tilt u. This issue is cleverly circumvented 
in [13] by instead considering the finite volume measure 



l*('\F(D*)°) = P 



.V\9D* 
D* 



almost surely. 
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on gradient fields on the torus, where v n is the uniform measure on the set of 
functions r\ : (Z$*)* — > R which can be expressed as the gradient of a function 
h: — y R. By construction, [i n is shift invariant, has tilt u and both of 
these properties are preserved in the limit as n — > oo. 

We will now explain how to use the method of dynamic coupling to prove 
that the gradient field of P^*, (p(b) = u • (yb — Xb) as before, converges to 
SEGGS n . Our proof will also yield an alternative construction of the Funaki- 
Spohn state in the special case u = 0. We will include a statement of this 
result here as well as a short sketch of the proof since this will be important 
for the proof of Theorem 1.1. We note in passing that Theorem 1.2 gives a 
much better coupling which will be critical for the proof of Theorem 1.1, but 
this result uses the convergence of the finite volume gradient fields hence we 
cannot simply apply it here. 

Fix a tilt u. Suppose that D n is any sequence of bounded domains in 
7i d converging locally to Z d . For each n, let r? n ~ PL, Suppose that 

7] ~ SEGGS U . Fix x n £ dD n and let h n ,h n ' S be the height fields associ- 
ated with the gradient fields r] n ,r], respectively, both set to vanish at x n . By 
the Brascamp-Lieb inequality (3.4), 

(4.6) Yar(h n ' s (x) - h n > s (y)) < Clog(l + \x - y\). 

Here, h n ' (x) — h n,s (y) is interpreted as Y^i=i v{h) where bi,...,b n is any 
sequence of bonds connecting x to y. Of course the same is also true with 
h n in place of h n ' S . Let R n = diam(L> n ). As h n (x) = Eh n ' S (x) for x E dD n , 
the Brascamp-Lieb (3.5) and Chebychev inequalities thus imply that 

(4.7) P[ max \h n ' S (x) - h n (x)\ > (logi? n ) 2 ] = 0{R- WO ). 

x&8D n 

Assume that (h™ ,S , hf) is the stationary coupling of h n and h n,s conditional 
on h n,s \g£) n . Then as /i" = h™' S — h™ satisfies the parabolic equation (4.3) 
and h is static on dD n , the maximum principle implies 

max |/i"(j;)| < max |/iq(x)|. 

x£D n xGdD n 

Combining this with (4.7) implies 

P[max \7%(x)\ > (logi? n ) 2 ] = O(R- 10 °). 

x&D n 

The Nash continuity estimate (Lemma B.2) applied to thus implies that 



(4.8) 



E[ max \Vh"(b)\] = O^R^ c ) 
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for e, C > fixed. This proves the desired convergence. 

Existence in the special case of u = can be proved in a very similar 
manner. The reason is that, in this case, h n ~ P^, hence Eh n (x) = 0. Thus 
it is clear that the subsequential limits of gradient fields, which exist by the 
Brascamp-Lieb inequalities, have zero tilt. The Brascamp-Lieb and Cheby- 
chev inequalities imply that the maximum of h n is with high probability 
0(logi? n ). Thus if h n denotes the law of the GL model on the domain D n 
given by shifting D n by one unit, then using the argument of the previous 
paragraph we can couple h n and h n such that V(h n — h n ) is with high prob- 
ability 0(R!n~^ NC ) at distance at least Rn, ( > 0, from both dD n and dD n . 
Therefore the subsequential limits of P^ n are shift-invariant which proves 
the existence of a zero-tilt shift-invariant Gibbs state. Uniqueness (and also 
existence of limits) follows by taking two such states rj, rj', then applying the 
argument of the previous paragraph. 

We have obtained: 

Theorem 4.4 (Equivalence of Ensembles) . 
// (D n ) is any sequence of bounded domains in Z d tending locally to Z d 
and, for each n, h n is an instance of the GL model on D n with boundary 
conditions ip n G B^-(D), then rf 1 = Vh n converges weakly to SEGGS U and we 
have 

E[ max \ V n (b)-r)(b)\] =0^(R € -^ NC ). 
beD«(Ri) 

4.3. Invariance under Reflections. The following proposition will be es- 
pecially useful for us later in the next section when applied to / = V". We 
let <f h ,ip v : R 2 —7- R 2 be the maps which reflect about the horizontal and 
vertical axes, respectively, and ^(b) = (p h {b — x), (f v x {b) = ip v (b — x) for 
x G Z 2 . 

Proposition 4.5. Fix a tilt u, x G Z 2 , and let f: R ->■ R be an even 
function. Let b h = tp%(b) for b G (Z 2 )*. If r, u ~ SEGGS U , then (/(r/ u (6)) : b G 

(Z 2 )*) = (f(ilu(b h )) : b G (Z 2 )*). The same is also true when the horizontal 
reflection is replaced with vertical reflection. 

PROOF. Let w = Lp h (u) and rj w ~ SEGGS W . Let s b = 1 if b is orien- 
tated vertically and —1 otherwise. Since / is even, we know that f(rj u (b)) = 

f(sbVu{b)). Since (si,r) u (b)) = (rj w (b)), we thus have that 
(f{riu(b))) = (f(r, w {b)))^f{ Vu {b h )). 
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The reason for the last inequality is that (r]u(b h ) ■ b G (Z 2 )*) is still a 
shift-ergodic Gibbs state but with tilt uu. □ 

5. Correlation Decay. Suppose that F C Z 2 is a bounded, connected 
domain with R = diam(F). Let £,£ G B^-(F) and assume that (hj,h^) is 

the stationary coupling of P^,,P^. Throughout, we let (3 = (5{u) as in the 
statement of Theorem 1.1. The main result of this section is that V"(V/i^(6)) 
and \7h(b), h = — , are uncorrelated when averaged against a A 13 - 
harmonic function. 

Theorem 5.1. Suppose that xo G F with dist(xo, dF) > R a+e for a, e > 
0, E = B(xq, R a ), and let g: F — > R be A@ -harmonic. We have that 

(5.1) E V"(Vh<(b))Vh{b)Vg{b) = ^ a u (b)B[Vh(b)]Vg{b)+ 

b£E* b£E* 

_ (jR£+a( l-p CD ) ||V5||oo) 

for p = Pcd(V) > and a u (b) = E[V"(r/(6))] with r] ~ SEGGS U . 

Note that a u depends on b only through its orientation (either vertical 
or horizontal) by the shift-invariance of n. Our typical choice of g will have 
|| V^Hoo = 0(R~ a ), in which case the exponent in the error term is actually 
negative. The idea of the proof is to show that replacing h(x) by its average 
h p (x) over the ball B(x,R p ) introduces a small amount of error. The ad- 
vantage of this replacement is that the time-derivative of h P t possesses more 
regularity than that of ht- This allows us to replace the left hand side of 
(5.1) with 

E V"{Vh c T {b))Vh (b)Vg(b). 

b&E* 

The proof is then completed by coupling /iy to jj ~ SEGGS U conditional on 
/io, which can be accomplished at the cost of negligible error by the argument 
used to prove Theorem 4.4. 

5.1. Change of Environment. Let rjt follow the SEGGS U dynamics inde- 
pendent from (/ij,/ij). That is, rjt solves the infinite dimensional SDS 



dVt((x, y))=[J2 V ^t(b)) - Y V'M&)) dt + V2(dW t (y) - dW t {x)) 

\b3y b^x I 
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for b G (Z 2 )* where Wt(x),x G Z 2 , is a family of iid standard Brownian 
motions; see Section 9 of [11] for a discussion of the existence and uniqueness 
of solutions to this equation. Fix T > 0, let ct(b) = V"(r] t (b)), and let p be the 
transition kernel of the random walk jumping with rates br-tip), t G [0, T], 
stopped on its first exit from F. 

Proposition 5.2. Suppose that we have the setup as Theorem 5.1. Let 
7i;72 £ (0,a] and5i = 4a— 4ji — jiPEC- Let S2 = R? 12 . There exists |i? 271 < 
Si < R 2 " 11 such that the following holds. Let p 1 ,]) 2 ~ p associated with 
independent environments rj 1 ,!] 2 ~ SEGGS U which are in turn independent of 
ho. Let S = Si + S 2 and let 

h t (x) = 'Yp 1 (S -t, Si; x,y)p 2 (0, S 2 ;y, z)ho(z), S - S 2 < t < S. 



y.z 



There exists a coupling of{r] ,rj ,ht) and (V/i, h) such that 

E (yh S (b) - Vh S (b)) 2 = 0_(^+2«+272-4 7 i+<52) + _(^+<5i) 5 
beE* 



E 



max sup \Vh^(b) — r]j( 

h ^ E * S/2<t<S 



O x (R e ~^ c ), 



for pec > depending only on V. 



The following heat kernel estimates are crucial ingredients in the proof of 
the proposition. 

Lemma 5.3. Suppose that x G F with dist(x , dF) > W +e for 7, e > 
and let E = B(x , K 7 )- For T = R 2 ^ and \T <t x <t 2 < T, we have 

Y \l{uM;x,y) - q(u,t 2 ;x,y)\ 2 = 0{\t x - t 2 \ 2 R^ 2 ^) 

x£E j/S-F 

for < u < T/4 and 
/■T/4 

\Vq(u,ti;b,y)-Vq( U ,t 2 ;b,y)\ 2 d U = 0(\ti-t 2 \ 2 tf- 2 ^c) 

where q is the transition kernel of a random walk on F evolving with rates 
av < dtib) < A v . 
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Proof. Using that dtq(u, t; x, y) = dt(b)Vq(u, t; x, 6), we have 
rT/4 



EE/ I *i; 6 ' f) - t2; 6 ' y)i 2du 

rt 2 rT/4 

<\h-t 2 \J2J2 / |3tVg(u,t;MI 2 *«# 



/•t 2 /-i/4 

<C|ti-t 2 | EE// |Wg(«,t;6, 
fees* 6'eF* ^° 



t 2 /-T/4 

?/)| 2 dudi, 



for C > depending only on V. Applying the Cacciopoli inequality (B.4) 
to the first time and spatial coordinates, we see that this is bounded from 
above by 

(5.2) (l)fo-ZM£ £ T [ T/2 \Vq(u,t;x,b')\ 2 dudt. 

The Nash-Aronson estimates (Lemma B.l) imply that the contribution 
to the sum given by those b' G F* with dist(6',F) > K y+e / 2 is negli- 
gible in comparison to the upper bound we seek to prove. For b' G F* 
with dist(6', F) < K 1+e / 2 , the Nash continuity and Nash-Aronson estimates 
(Lemmas B.2, B.l) imply 



\Vq(u,t;x,b')\ = 0{R 



for < u < T/2 and t\ < t < t 2 . Inserting this into (5.2), we arrive at the 
bound 

\h - t 2 \ 2 0(R^ +e ) ■ 0(iT 27 « NC - 47 ) = 0(\h - t 2 \ 2 R e ~ 2 ^ c ). 

This proves the second part of the lemma. The first is exactly the same 
except the application of the Cacciopoli inequality is unnecessary. □ 

Lemma 5.4. Suppose that x G F with dist(x , dF) > R? +e for 7, e > 
and let E = B{x ,R?), E = £(x , i? 7+e/2 ), E 1 = B(x ,R~> +e ). Let q,q' 
be the transition kernels for two random walks on F jumping with rates 
&v < dt(b), d' t (b) < A\>, respectively. With T = R 21 , assume that dt = d' t for 
all <t < gT and let 

doo = sup max \d t (b) — d' t (b)\. 

0<t<T b ^ E t 
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Uniformly in x £ Eq we have that 

Y, \q(0,T;x,y)\ 2 +Y F |Vg(0, t; x, b)\ 2 dt = O^BT^) 

y£E{ b&E* J ° 

where q = q — q' . 

PROOF. By definition, 

d t q(0,t;x,y) = [£ t q(0, t; x, -)](y) and <9 t </(0, t; x, y) = [C' t q'(0,t;x,-)](y) 
where 

Ctf(y) =£dt(i)V/(6) and C'J{y) = ^4(6)V/(6). 
Consequently, 

(5.3) d t q 2 = 2q(C t q + C t q') 

where Ct = Ct — C' t . Using the same proof as the energy inequality (Lemma 
4.1), by integrating both sides of (5.3) from to T then using summation 
by parts, 



Y \mT;x,y)\ 2 + 2 Y [ d t (b)\Vq(0,f,x,b)\ 2 dt 



yeE! b£E* 
i-T 



< 



bedEl Jo 
(-T 



d t (b)\q(0, t; x, x b )Vq(0, t; x,b)\ + \d t (b)q{0, t; x, x b )Vq'(0, t; x, b)\ 



dt+ 



2 V f \d t {b)Vq(0,t;x,b)Vq'(0,t;x,b)\dt. 
best Jo 

By the Nash-Aronson estimates (Lemma B.l), 

[ \q(0,t;x,x b )Vq(0,t;x,b)\dt = O(R 3 ^ +e exp(-c'R € )) = 0(exp(-cR £ )) 

b£dE* ^° 

for some c, d > depending only on ay , Ay ■ The other boundary term is 
similarly of order 0(exp(— cR e )). Consequently, 

Y \q(0,T;x,y)\ 2 + £ F | Vg(0, t; x, b)\ 2 dt < 0(exp(- Cj R e )) + 
yEEi beE{ Jt I 2 

(5.4) Cdoo V f \Vq(0,t;x,b)Vq'(0,t;x,b)\dt. 
beEl Jt / 2 
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The reason that the lower bound of integration is T/2 rather than is dt = d' t 
for all t<T/2. If 

f \Vq{0,t;x,b)\ 2 dt = 0(exp(-cR e/2 )) 

b£E* ^ T / 2 

then we are obviously done. If not, we apply Cauchy-Schwarz to the sum on 
the right hand side of (5.4) and rearrange to arrive at 

f \Vmt;x,b)\ 2 dt<CdloYl f \Vq'(0,t;x,b)\ 2 dt, 
b£E* •' T / 2 beE* ^ T / 2 

increasing C if necessary. By the Nash-Aronson estimates, 

\q'(0,T;x,y)\ 2 = O(R- 2 ~<) and 

ye Ei 

rT 

Y / \q'(0,t;x,x b )\\Vq'(0,t;x,b)\dt = 0(exp(-ciT)). 
bedEl Jt I 2 

Therefore, by the energy inequality (Lemma 4.1), 

(5.5) V [ T \Vq'(0,t;x,b)\ 2 dt = O(R-^), 

r-z. Jt/2 



which in turn implies 

(5.6) V f T \Vq{Q,t-xM 2 dt = 0(d 2 oa R- 2 "<). 

bEE* Jt / 2 

Inserting (5.5), (5.6) into (5.4) proves the lemma. □ 

Lemma 5.5. Suppose that we have the same setup as Theorem 5.1 and let 
p, r\ be as in the introduction of this subsection. Let S = R 2j for < 7 < a 
and 

h t (y) = YP(S-t,S;y,z)h (z) forO<t<S. 
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There exists a coupling of (rj,h) and (V/i^/i) such that 
r s 

(5.7) ET \Vh t (b)-Vh t (b)\ 2 dt = O x (R e+4a - 2 ^ pEC ), 

C~±, J3S/4 



b&E 



(5.8) E Y \hs(x) - h s (x)\ 2 = O x (ir+^2 7 - 7PEc; 



(5.9) E 



max sup \Vh^(b) — r] t (b)\ 

h&E * S/2<t<S 



where p^c > depends only on V and r] is independent of (V/iq, ho). 

This constant pec from Proposition 5.2 is the same as that appearing in 
the statement of this lemma. 

Proof. Let T = R 2a and let 770 ~ SEGGS„ independent of (h^ho). As- 
sume further that the evolution of the Brownian motions driving rjt in F are 

independent from those of (h^,h^) until time T — S, S = R 2 " 1 , after which 
they are the same. Let ct(b) = V"(rjt(b)), ct(b) be as in (4.2) with h^,h^ in 
place of h^,h^, and let p be the transition kernel of a random walk in F 
stopped on its first exit jumping with rates CT-t(b). Note that 

h t (x) = ^p(T-t,r ;3 ;,y)fro(y) for < t < T. 

y&F 

Set S' = S — R 2a , a G (0,7) to be determined later, and define environments 

(c t (b) for T-S<t<T ^ fct(6) foiT-S'<t<T 

c# = < , cAb) = < 

\c t {b) for <t <T - S \ct(b) for < t < T- S' . 

Let p, p' be the transition kernels for the random walks in F stopped on 
their first exit jumping with rates CT-t(b) ,c^ T _ t (b) , respectively. Finally, let 
p = p—p and p' = p—fi. For S < t < T, we have p(S, t; x, y) = p(S, t; x, y) = 
p 1 (S, t;x,y), hence 

Y ^(Vp( U ,r;6,2/)) 2 = Y, E (E^^MP^T;^^ . 

b&E* yeF beE* y£F \z&F ) 

By Jensen's inequality, this is bounded from above by 
4 E E ((yp(u, S; b, z) - Vp(u, S'; b, z)) 2 

b£E* z£F ^ 

+ (Vp'(n, S'; b, z) - Vp'(u, S; b, z)) 2 + (Vp'(u, T; b, z)) 2 ^j = h + I 2 + I 3 . 
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Fix a base point ao € dF, set hf(ao) = 0, and let /if solve V/if (6) = ryt(fe). 
Applying the Nash continuity and Nash-Aronson estimate (Lemma B.2) to 

/if — and h\ — h^, similar to the proof of Theorem 4.4, yields 

(5.10) doo = sup sup \cr-t{b) ~ c T -t(b)\ = O{M IT a ^c) 

0<t<S' b£E* 

where M t = \\ht\\oo + \\ht-h?\\ oo, as V" is Lipschitz. Here, we are taking the 
maximum over x S F. Let q(s, t; x, y) = p(T—t, T—s; y, x) and q'(s, t; x, y) = 
jf(T — t,T — s;y,x). Then q, q' are the transition kernels for random walks 
jumping with rates q(6),c^(6), respectively. The previous lemma thus yields 

E E [ T (W(n,T;b,y)) 2 du = 0(M^R^ 2 - 2 ^ ) 

b£E* y£F ^° 

since the contribution to the sum given by those y £ F with dist(y,.E) > 
j^a+e/i - g negligible. Since we can cover E by 0(R 2( - a ~^) disks of radius W , 
applying Lemma 5.3 to the terms corresponding to I±, I2 gives us the bound 

EE/ (Vp(u,T;b,y)) 2 du 

(5.11) =0 ( jR ,/2+4 ( 7+2( a - 7 )-2 7 € N c) + 0{M$R e / 2 - 2a ^ c ). 

Observe that M is the only random quantity on the right hand side. Let 
h t (x) = J2p(T-t,T;x,y)h (y) for < t < T. 

y&F 

Note that 

ht{x)= E P( T ~ ^^jM^^y.^oW 

y,zeF 

= Y J P(T- t, S; x, y)h T ~s(z) for T - S < t < T. 

yeF 

Hence as t 1— > p(T — t, S; x,y), T — S < t < T, is independent from hx-s, 

it follows that /H+er-S) = ht, < t < S, with h as in the statement of the 
proposition. We can write 

EV f T (Vht(b)-Vh t (b)) 2 dt 
b&E * Jt-s/a 

=E E f T \TVp(T-t,T;b,y)h T _ t (y)) dt. 
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The terms in the summation over y £ F which are of distance at least 
f rom £ ma k e a negligible contribution to the summation by the 
Nash-Aronson estimate. Consequently, by making a change of variables and 
applying the Cauchy-Schwarz inequality, we see that it suffices to control 



(5.12) E 



r S/4 

R 2a+e J2"£ (Vp(«, T; b, y)) 2 M 2 du 

Ucl TP* 77 J 



feS-B* y£F 



Choosing a = 7£nc/4 and applying (5.11) yields that the expression in 
(5.12) is bounded by 

sup ( EO(M*R 3 ^ 2+4a -^-^) + EO(M$M*R 3 ^ 2+2a - 2 ^c/*)) . 

0<u<S/4 V / 



Now, Lemma 3.3 and the Cauchy-Schwarz inequality imply E(M 2 +Mq = 
0-^{R e l 2 ) uniformly in u. Thus, we are left with the bound 

_^2 e +4a-2 7 - 7 ^ c /2^ 

Taking p^c = £nc/2 gi yes (5-7). Equation (5.8) follows from exactly the 
same argument except using the first part of Lemma 5.3 rather than the 
second. The final part of the proposition is immediate from the construction 
and the Nash continuity estimate. □ 

We can now prove Proposition 5.2. 

Proof of Proposition 5.2. We now construct couplings as follows. 
First, we couple (rj 2 ,h 2 ), (Vh, h) as in Proposition 5.2 for 7 = 72. Equa- 
tion (5.8) implies that with S2 = R 2 ^ 2 we have 

(5.13) E £ \h S2 (x) - h 2 2 (x)\ 2 = O^R^ 2 ^ 2 ). 

where E = B(xo,2R a+e / 100 ). Now we apply Proposition 5.2 a second time 
except with 7 = 71 and starting at S2 to yield a coupling ((77 , h ), (Vh, h)). 
Equation (5.7) implies the existence of |i? 271 < S\ < R 2 " 11 such that with 
S = Si + S2 we have 

(5.14) E £ |VM&) - V^(6)| 2 = G\(iT+ 51 ). 

b£E* 
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Let p l be the kernel associated with rf . In this coupling, ht for < t < S2 
is independent from p 1 and p l ,p 2 are independent. Let 

h(x) = 52 P l {S -t,Sr,x 1 y)p 2 (Q,S2)y,z)hQ(z) 

y,z£F 

= ^P 1 ^ - t, Si;x,y)hs 2 (y) 

ydF 

for S 2 < t < S. We have 
EJ2 \Vh S (b)-Vh l s (b)\* = Ej2 i^P^'^yWsM -hsM) 

beE* b£E* \yeF 

^ I E E^^ ' 5 !!^)) 2 ( E E (^ 2 ^)-^ 2 (y)) 2 ) +0(exp(- J R 

where the last inequality came from Cauchy-Schwarz and the Nash-Aronson 
estimates (Lemma B.l). It follows from equation (1.4) of Theorem 1.1 from 
[6] and the Nash-Aronson estimates that 

52 52 V(VpH0,S i; b,y)) 2 = 0(R 2a+ ^). 

Combining everything proves the proposition. □ 
5.2. Approximation by the Average. 

Proposition 5.6. Suppose that we have the same setup as Theorem 5.1. 
There exists pa > depending only on V such that the following holds. If 
h p (x) is the average of h{x) on the ball B(x,R p ), then 

E 52 V"(V4 2 (6))(V%(6) - Vh p Tl (b))Vg(b) 

b£E* 

(5.15) =OA(R e+p+a(1 - pA) \\Vg\U 
for T\ <T 2 . 

Proof. While the proof for T% ^ T 2 does not introduce any additional 
technical challenges, in the interest of keeping the notation simple we will 
only provide the proof for the case T± = T 2 = S, with S from Proposition 
5.2. Let 

( l+^PEc\ , 

71 = a f an( l 72 = a - 

\ 1 + JPEG 



10 
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Let 61,62 be as in Proposition 5.2 with these choices of 71,72- Note that 
then there exists pa > depending only on V such that 

61 < —apA and 2a — 471 + 62 + 272 < —apA- 

Consequently, by Proposition 5.2 it suffices to prove (5.15) with h,h p in 
place of h,h P where h P [x) is the average of ht{x) on B(x,R p ). Moreover, 
(5.9) combined with Lemma 4.3 imply that it suffices to prove (5.15) with 
V"(v4(k)) replaced with c l {b) = V"(r/^). For /: F -»• R let 

Cf(x) = Y J c\b)Vf{b). 

b9x 

Summation by parts implies it suffices to bound 

(5.16) &1 (b)(hs(x b ) - h p s (x b ))Vg(b) - £>s(*) " h p (x)]Cg(x). 

bedE* x£E 

By the Nash continuity estimate (Lemma B.2), we know that 
V\hs(x b ) - h p s (x b )\ = O x (R e+ ( p - a ^ c ), 

hence the boundary term in (5.16) is of order 0^(R e+ ^ p ~ a ^ NC+a \\ V<?||oo) = 
0^(i? e+a ^ 1_PA ^ || V5H00), shrinking pa if necessary. 

We now deal with the interior term. For y,6 € Z 2 , let y e = y + 9. We 
are going to omit the times when referring to p 1 and just write p l (x,y) for 
p l (0, Sf,x,y). Say that a bond b is "positively oriented" if it points either 
in the positive hortzonal or vertical directions. For each triple (x,y,b), let 
y b = y if b is positively oriented and 2x — y otherwise. The latter is the 
reflection of y about x. With B p = B(0,R p ), we can rewrite the interior 
term of (5.16) as 

' p ' x<=E b3x y,z£F 8eB p 

Using the independence properties of rf,ho as well as the A^-harmonicity 
of g, we see that this is the same as 

j^-j Yl w(x,y,b,9)B[p 2 (y b ,z)-p\x,z)]B[h (z)]^g(b) 

p x,b,y,z,6 
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where the summation is over x£E,b3x,y,z£F,9£ B p and 

w[x, y, b, 9) = E[c\b){p\x, y b ) - p\x e b , y b ))] . 

Here we are crucially using that p 2 (x, z) does not depend on y. For b 3 x let 
b' 3 x have the opposite orientation of b. By Proposition 4.5, we have 

w(x,y,b,9) = w(x,y,b',9) 

Consequently, we can rewrite our sum as 

w(x,y,b,9)rE[{p 2 (y b ,z)-p 2 (x,z)) + (p 2 (y b ,,z)-p 2 {x,z))]Vg(b) 

" x,b,y,z,9 

+ E[p 2 (y b ,, z) - p 2 (x, z)} (Vg(b') - Vg(b))\ E[h (z)} 

=W~\ S y ' 6 ' ^ ( A ( x > y > z > b ^ 9 ^ + B ( x > y > z > b )( V 9(b') - Vg(b)) J EpioCs)] 

' p ' x,b,y,z,6 ^ ' 

where the summation is now only over positively oriented bonds. 

We will deal with the term involving A first. By the Nash-Aronson esti- 
mates (Lemma B.l), the sum over y of distance from x more than W 1+e 
from x is negligible. Similarly, we may ignore those z with \z — x\ > W 2+e . 
For \y — x\ < i2 7l+e , it is a consequence of Theorem 1.1 equation (1.5b) of 
[6] and the Nash-Aronson estimates (Lemma B.l) that 

(5.17) \A(x,y,z,b)\ < CR 2 ^ +2€ ~^ 2 . 

Indeed, this can be seen by rewriting the difference as a sum of 0(R 2ll+2e ) 
discrete second derivatives. Using that Eho(z) = 0-^(R e ) from Lemma 3.3, 
we thus have 

j— ^ w ( x ,yib,9)A(x,y,z,b)Vg(b) 

" x,b,y,z,9 

=nr\ 2 ^y^o^R^^nvgU). 

P x,b,y,6 

Now, equation (1.4) of Theorem 1.1 of [G] implies 

Ej2H^y,b,0)\ = o(Rp-^) 

y 

uniformly in x since \x — x b \ < R p . Since the sum over x includes 0(i? 272 ) 
terms, our total error is O^(R 10e+p+l1 || Vg||oo). 
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We now turn to the term involving B(x, y, z, b). Since g is A^-harmonic so 
isx^ Vg((x,x + ei)) hence we have that Vg(b')-Vg(b) = 0(R^ 2 \\V g\\oo). 
Again applying equation (1.4) of Theorem 1.1 of [6] we thus see 

B(x,y,z,b)(Vg(b') - Vg(b))E[h (z)] = O^R^-^WVgU). 

This is of an even smaller magnitude than the corresponding term with A, so 
we also get an error of O^-(i? 10e+p+71 1| V<?||oo). The proposition now follows 
from the explicit form of 71 . □ 

5.3. Change of Time. 

Proposition 5.7. Suppose that we have the same setup as Theorem 5.1. 
Let T = R 2 ' y forO<-f<a and fix g: F — > R. We have that 

(5.18) E ^ V"(Vh C T (b))(Vh T (b) - Vh (b))Vg(b) 

b&E* 

=O x (ii e+Q ( 1 -^oT)+7||v 5 || 00 ) 

where pcoT > depends only on V. 

Proof. For p > 0, Proposition 5.6 implies that replacing hx,ho by 
h^, /iq, respectively, in (5.18) introduces an error of 

(5.19) Ox( J R e+p+(1 - M)a ||V 5 || 00 ). 
We will now prove that 

E Y, V"(v4(6))(V^(fc) -Vh p (b))Vg(b) 

b&E* 

(5.20) =o x (ri? e+a -^|v 5 || 00 ). 

Since Ylb&E* (^9(b)) 2 = 0(R 2a \\Vg\\ 2 XJ ), applying the Cauchy-Schwarz in- 
equality to the expression on the left hand side implies that it suffices to 
show 

E £ (Vh P T (b) - Vh P (b)) 2 = 0^(T 2 R^). 

b&E* 

As 

\d t h P t {x)\ < H £ | V ^(6)| 

bedB*(x,RP) 
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we have that 



E (V/4(&) - V/io(6)) 2 < CiTE f Y \Vdth[ 

b£E* ^° b£E* 

j-i* I,^ 77* 



where £7i = B(xq, R a + R P ) and the final equality comes by stationarity. We 
know that the latter quantity is of order O^R 6 ) by Lemma 3.3. Equating 
the exponents in (5.19), (5.20) gives the equation 

p + (1 - p A )a = 27 + a - p, 

which leads to the choice 

p = 7 + — a. 

Combining everything gives the proposition, where pcdT = 5 PA- D 

5.4. Proof of Theorem 5.1. Assume that T = R 2 ^ where 7 = apcoT/2- 
By Proposition 5.7 it suffices to estimate 

(5.21) Yl E[V // (V4(6))V4(6))]V 5 (6) 

fees* 

provided we pay an error of 0^;(i? e+a ^ 1_PCoT ^ 2 ^ || V^Hoo). Now the result fol- 
lows from an argument similar to the proof of (5.10) from the proof of 
Proposition 5.2. Indeed, we let rjt follow the SEGGS n dynamics driven by the 

same Brownian motions as (h^,h^) but with 770 independent from (/iq,/iq), 
then use the Nash continuity estimate (Lemma B.2) to argue that 



P[max|V/4(6) - rjr{b)\ > fi £ " 7&c ] = O x (iT 



100 \ 



Thus applying the Cauchy-Schwarz inequality, we see that replacing V"{Vh C T (b)) 
with V"(rjT(b)) in (5.21) introduces an error of order 

n(Vho(b)) 2 }) • O x (iT||V<?Hoo) • O x (i?^ NC ). 

\b£E* J 

The result now follows as the first term is of order 0^(R e ) by Lemma 4.3. □ 
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6. Harmonic Coupling. Throughout this section we will make use of 
the following notation. For F C Z 2 and (f>: dF — > R, we let be the 
law of the GL model on F with boundary condition eft. We will denote by 
a random variable distributed according to P^, where F is understood 
through the domain of definition of (f). Finally, for g : F — > R we let Qp 9 be 
the law of (h^—g). We will write E^ for the expectation under P^ if we want 
to emphasize ip and, similarly, E^'^ for the expectation under a coupling of 

P^, Pj£ if we want to emphasize both if) and ip. We also let ft = (5{u) as in 
the statement of Theorem 1.1 for a fixed tilt u 6 R 2 . 

Suppose that /i, v are measures with \i absolutely continuous with respect 
to v. Recall that the relative entropy of \i with respect to v is the quantity 



H( M |i/)=E„ 



log 
du 



We begin by fixing D C Z 2 with R = diam(D) < oo. Fix A > and let 
ip,ip £ B^-(D). Morally, the idea of our proof is to get an explicit upper 
bound on the rate of decay of the symmetrized relative entropy 

H(P^|Q^) + H(Q^|P^), 

where h is the A^-harmonic extension of ip — tp from dD — > D, as R — > 
oo, then invoke Pinsker's inequality, the well-known bound that the total 
variation distance of measures is bounded from above by the square-root of 
their relative entropy [5]: 

(6.1) || M -z/|| 2 y < iH( M |i/). 

We will show shortly that the symmetrized relative entropy takes the form 

(6.2) Yj ^c(b)Vh(b)(Vh(b) - Vh(b)) 

beD* 

where h = W 3 — h^ and c(b) is a collection of conductances which are random 
but uniformly bounded from above and below. In the Gaussian case, c(b) = c 
is constant, hence one can sum by parts, then use the harmonicity of h to 
get that the entropy vanishes. The idea of our proof is to use Theorem 5.1 
repeatedly to show that this approximately holds in expectation: 



Y B^c(b)Vh(b)(Vh(b) - Vh(b)) 

b&D* 

(6.3) = a u (b)E^Vh{b)(Vh(b) - Vh(b)) + O^iT" 5 ) 



beD* 
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for some 5 > 0, where a u (b) = E[V"(t?(6))] for rj ~ SEGGS U . Note that a u (b) 
depends only on V, u, and the orientation of b. 

Theorem 5.1 is only applicable if the distance of b to dD is ^l(R^). This 
will force us to deal with a boundary term, the magnitude of which will in 
turn depend on the regularity of both V/t and h near dD. Since we make 
no hypotheses on ip, ip other than being pointwise bounded it may very well 
be that neither h nor h possess any regularity near dD. 

We will resolve this issue by invoking the length-area comparison tech- 
nique of Section 4. This gives us that, for e > fixed, there exists R 1 ~ e < 

R X D < R 2 D = Rd + Rl ~ 5e ^ 2Rl ~ e such that 

£ \Vh(b)\ 2 = O^R-*) 

where D(R 1 ,R 2 ) = {x G D : i?i < dist(x,dD) < R 2 }. Let g be the A 13 - 
harmonic extension of h from dD(R l D ,R 2 D ) to D(R l D ,R 2 D ). Note that g is 
the minimizer of the variational problem 

g^ Yl a u (b)(Vg(b)) 2 , g(x) = h(x) for x G 8D(R\), Rjj). 

beD*(Ri D ,n 2 D ) 

Indeed, the first order conditions for optimality are exactly that A^g(x) = 
for x G D(R l D ,R 2 D ). Consequently, 

E** £ |V 5 (6)| 2 = O x (i?-). 

Going back to (6.2), by invoking Pinsker's inequality, this implies that we 
can construct our initial coupling so that h is harmonic in D(R l D ,R 2 D ) on 
an event U with P[H] = 1 - O^iT^ 2 ). On we have that Vh(b) = 
O^R 66 ' 1 ) uniformly in b G dD(R D ) where R D = R l D + \R l ~ bt . Thus with 
high probability h has plenty of regularity a bit away from dD while ip — ip 
need not have any. 

Moving to the subdomain D(Rd) from D is also useful since it possesses 
the r-exterior ball property for r = Rq. This means that for every x G 
dD(R D ) there exists y <£ D(R D ) such that B(y,R D ) n D(R D ) = and 
x G dB(y, Rd)- The importance of this property is that it implies pointwise 
regularity of harmonic functions in D(Re>) near dD(R£>), more so than one 
has for such functions in D near dD without further hypotheses. This is 
related to the notion of "stochastic regularity" [16] and that random walk 
"exits much more quickly" from such domains. 
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(a) Domain without the ex- (b) Domain with the r- 
terior ball property exterior ball property 



Fig 2. The domain on the left hand side does not have the r -exterior ball property for 
due to the fjord in its upper right corner. However, the inner domain D(r) — {x £ D : 
dist(rE, dD) > r} shaded in light grey on the right hand side trivially does possess this 
property. The distinction is important since the r-exterior ball property is related to the 
regularity near the boundary of discrete harmonic functions. 



The rest of this section is organized as follows. In subsection 6.1, we will 
justify (6.2). The purpose of subsection 6.2 is to prove a form of (6.3). 
Finally, in subsection 6.3 we will put everything together to prove Theorems 
1.2 and 1.3. 

Before we proceed, we would like to emphasize that in this section we 
will prove that a certain symmetrized relative entropy decays like a small, 
negative power of R = diam(Z)). Because of this, many of our estimates will 
be accurate only up to very small powers of R. In particular, we do not try 
to derive the "best possible" exponents and make many cautious choices in 
order to avoid carrying around overly complicated exponents. 

6.1. The Symmetrized Relative Entropy. 

Lemma 6.1. Suppose that F C Z 2 is bounded and £: dF — > R are 
given boundary conditions. If g: F — >• R is any function such that g\gp = 
C — C then 

h(pS|q£ 9 ) + h(q£ 9 |pS) 

V"(Vh<(b))Vg(b)V(9 ~ h)(b) + 0((\Vh(b)\ 2 + \Vg(b)\ 2 )\Vg(b)\) 



beF* 

where denotes the expectation under any coupling o/Pi,,Pl and h 
hi - h<. 
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Proof. The densities p,q = q g of and Qp with respect to Lebesgue 
measure are given by 

p(/ l ) = ^exp(-^V(V(/ i VC)W)), 

p V b<=F* J 

q(h) = ^exp ( - V(V[(h + g)V(](b))) . 

With ~Er-< 9 the expectation under Q^ 9 , we have 

V((VhV ()(b)) -V(V((h + g)V ()(b)) 

V Jo / 



H(Q^|P^) + log|^= ^E<- 



6GF* 



Similarly, with E^ the expectation under Pi, 



H(Pj.|Q^)+log^= J>< 



q beF* 



V(V((h + g)\/ 0(b)) -V((Vh\/ 0(b)) 



, cP * V Jo / 



b€F 

Thus 

H(PJ|Q^) + H(Q^|P^) 

= V (V / V'(V(/i + s 5 )(6))da - E< / V'(V[(/i + (a - l)s)](&))<fc) V 5 (6) 
beF » V Jo Jo / 

6eF * V Jo Jo / 
Vg(b)V(g-h)(b). 

As V" is Lipschitz, 

f [ V'tV^ + (s-l)g)(b) + rV(g-h)(b))drds 
Jo Jo 

=V"(V/t c (6)) + 0(Vh(b)) + 0(Vg(b)). 
The lemma now follows by an application of Cauchy-Schwarz. □ 
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6.2. Estimating the Entropy. Suppose that ECZ 2 with diam(i?) = R 
and £, £ 6 B^-(E'). We are now going to give a general estimate of the 
symmetrized relative entropy of the previous lemma when g : E — > R, is the 
A^-harmonic extension of C = C — C horn dE to E. The error is going to 
be a function of the regularity of g, £, and the number of balls required to 
cover annuli near dE. When all of the boundary data is smooth, the error 
is actually negligible. To this end, we let 



If 



max 



IC(x)-C(v)l 



x,yedE \x — y\ 



Fix e > 0, let 7^ = e/c, M be the largest integer so that 7m < 1, and be 
the number of balls of radius R lk necessary to cover E(R" lk , R" /k+1 ) = {x G 
E:R?*< dist(x, dE) < Finally, with 1 < I < M fixed, let 

£ E = \E*\\\Vg\\l +£\\Vg\\ 00 , 
4 = l(^T(^)ll|V<7||oo^, 

M 

e\= ^ e+7fc(1 - pcD) llv 5 ||oo, 

k=e+i 



where 



£ e = (i? 1/2 ||Cllf + ^- 1/2)PB ||CIU) and £ = IdE*]^^ 1 



Proposition 6.2. We have that 

H(P< |P<) + H(P| \P%) = Oj(£ E + £ I + £ B ). 

Proof. Let (h 1 *, h£) be the stationary coupling of P E , P E and h = h£ — 
h^. We begin with an a priori estimate on the Dirichlet energy of h. First 
fix b g dE*. Lemma B.4 implies that Vh(b) = OiR^WCWv + ir^/ 2 ||C||oo)- 

Applying Lemma 4.1 to the stationary dynamics (h^, hf), we have that 
£ |VM&)| 2 <G £ |C(a*)||VE(6)|. 

b£-E* feeds* 

Consequently, we have 



(6.4) 



fees* 
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We now break the right hand side in the statement of Lemma 6.1 into 
three terms: 



(6.5) ^E«(|V 5 (6)| 2 + |VM&)| 2 )|V 5 (6)|, 

V"(VhHb))Vg(b)V(g-h)(b) 



b&E* 



(6.6) Y. E ° 

be{E')*{R?i) 

(6.7) Y, E ° 

b£E*(FOi) 



V"(Vh<(b))Vg(b)V(g-h)(b) 



Estimate of (6.5) 

The first term in the summation is trivially bounded by |£7*|[|V<7||^. By 
(6.4) the second is at most ||V<?||oo£- This gives an error of 0(£e)- 
Estimate of (6.6) 

The error is easily seen to be | (E')*(R ye )\ \\Vg\\oo 5 where 5 = max xg£ ;/( K 7£) \g{x) — 
h(x)\, so we just need to estimate 5. Fix x € E'(R? e ). We know that we can 
write g(x) = E£(X T ) where X is a random walk initialized at x with bounded 
rates and r its first exit from E. By Lemma B.4, the probability that X T 
exits at distance less than R 1 / 2 from x G E'(W l ) is 1 - 0(i? (7 ^ 1/2)pB ). 
Now, h admits a similar representation though the random walk has time- 
varying rates. Nevertheless, the same statement still holds. Consequently, 
\h(x) — g{x)\ < £ e , which leads to the desired bound. This gives an error of 

o(4). 

Estimate of (6.7) 

We break the summation over E* (R? e ) into the annuli E^ = E* (i? 7fc , R lk+1 ) . 
Each annulus can be covered by balls of radius R lk by hypothesis. Let B 
be such a ball. Applying Theorem 5.1 on B yields an error of 0^(i? e+7fe ( 1 ~' 9 cD) || Vg|| c 
Repeating this on all iVjt balls for t + 1 < k < M implies that (6.7) is equal 
to 

M 

Y a u (b)B[V(h - g)(b)}Vg(b) + Y N^R^ 1 ^ \\Vg\U 

b&E*{R~<e) k=l 

where a u (b) = E[V"(r/(6))], r\ ~ SEGGS M . Note that the error term is precisely 
£j. By our bound of (6.6), we can rewrite the above expression as 

Y a u (b)B[V(h - g)(b)]Vg(b) + 0^£ E + £j). 

b&E* 
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(a) Stage 1: Langevin Dy- (b) Stage 2: Harmonic in 
namics in D D(Rj D ,Rj D ) 

Fig 3. The entropy estimate consists of two stages of coupling, indicated by the images 
above. The annulus surrounded by a dashed line is D{R 1 D , R 2 D ), and inner light grey region 
in (c) is E = D(R D ). 

By summation by parts, we see that this is equal to 

Ep - g)(x)]^g{x) + 0^(£ E + £j) = C^(£ E + £j), 

b&E* 

where (3 = fi{u) as g is A^-harmonic. This gives an error of 0(£f). □ 
6.3. Proof of Theorems 1.2 and 1.3. 

6.3.1. The Initial Coupling. In this subsection we are going to show that 

we can construct a coupling of P^, , P^> so that with high probability the 
error terms from Proposition 6.2 are with high probability negligible when 

applied to P^,P^ where E = D(Rd), Rd = cR 1 ~ e ° , some €d > 0, and 

((, C) = (h^, /i^)|(aE)2. We will accomplish this using the following steps: 

1. Take the stationary coupling (h^, h^) of P^,, ~P D - 

2. Invoke Lemma 4.3 to find an annulus D(R D , R D ) on which the Dirich- 
let energy of h is controlled. 

3. Use Lemma 6.1 to show that we can recouple the laws on D(R D , R D ) so 
that with high probability h is A^-harmonic. On this event we will have 
all of the regularity that we need on dD(Rry), where R D < Rd < R?d- 

Fix en > so small that €d < (10~ 100 pcd A pe A £nc) 2 > where £nc 
is the exponent from the Nash continuity estimate (Lemma B.2), p-g is 
the exponent from Lemma B.4, and /?cd is from Theorem 5.1. We assume 




(c) The domain E — 
D(R D ) 
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R 1 ^ <R l D < 2R 1 ~ 6D has been chosen such that with R 2 D = R 1 D + R l ~^ D 
and R D = R 1 D + ifi 1 " 5 ^: 

(6.8) Yl |VM&)| 2 = O x (iT^), 

beD*{B} D ,Rl) 

(6.9) ^ \Vh(b)\ 2 = 7[ (R 3€D ). 

beD*(R D ) 

That such a choice is possible is ensured by Lemma 4.3. Let E = D(Rd). 

In order to apply Proposition 6.2 we need to make sure that we can 
arrange for the number of balls required to cover annuli near the boundary 
is not too large. Such estimates would come for free if D was a lattice 
approximation of a smooth domain in R 2 . As we want Theorem 1.2 to hold 
for general bounded subsets of Z 2 , such estimates do not necessarily hold 
uniformly but only on the average provided we are far enough from 3D. Let 
7/t = keo and let be as in Proposition 6.2. In particular, by using the 
averaging technique of Lemma 4.3 we can arrange for Rjj to be such that 

(6.10) N k = 0(i? 1+72 - 7fe ), \E(R yk )\ = 0(R 1+ ~ fk + 2 ), \8E\ =0(R 1+ ~ tl ) 
for all k < M, M the largest integer such that jm < 1. 

Lemma 6.3 (Harmonic Coupling at the Boundary). There exists a cou- 
pling (h^,h^) o/pJ,p£ such that 

U = (h = - ift is A? -harmonic in D(R]j, R 2 D )} 
occurs with probability 1 — 0^(R~ tD / 2 ). 

Proof. Let h = h^ — h$, F = D{R} D ,R 2 D ), and let g: F -> R be A 13 - 
harmonic in F with boundary values h on OF. By our choice of Rd, 

^E^(Vg(b)) 2 = O^R-^) 

b&F* 

as g is harmonic in F, has the same boundary values as h, and h satisfies 
the same estimate. Let C, C = (^,^)| OFxdF- Conditional on (C,C)> let 
Pp, P F have the laws of the GL model on F with boundary conditions £, £, 
respectively, and let Qp 9 have the law of h** — g where h£ ~ P^. It follows 
from the Cauchy-Schwarz inequality and Lemma 6.1 that 

E^[H(P F |Q£ 9 ) + h(q£ 5 |pJ)] = o x (R-*°). 

The lemma follows by invoking Pinsker's inequality (6.1). □ 
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6.3.2. Regularity Estimate. In the following lemma we will use (h^, h^) 

to indicate a random variable with joint law given by the coupling of P^, P^ 

from Lemma 6.3 and H the corresponding event. Let (£, £) = (/i , h^)\gExdE- 
Let g : E — > R be the A^-harmonic extension of £ = £ — £ from <9i? to -E7. 
Recall the definition of ||C|lv J us ^ before the statement of Proposition 6.2. 

Lemma 6.4. There exists 1 < cp < 10 so that 

(6.11) (E^dKllf nn)^ = O^R^-- 1 ), 

(6.12) (E^^[max|V 5 (6)|n w ]y /P = O^R^" 1 ) 
for every p > 1 . 

Proof. By construction, E has the r-exterior ball property for r = 
ft}— zd Furthermore, if x,y 6 with |x — y| < |i? 1_5eD then the shortest 
path connecting x to y in Z 2 is contained in B(x, |i? 1_5eD ). Consequently, 

(6.13) |C(x)-C(y)| <g\x-y\ 
where 

# = max{|Vh(6)| :b€D*(R 1 D + \R l ^ e ° , R\> + f R 1 ^ )}. 

Let M = max{|/i(x)| : x S £>}. Since h is harmonic in D(Rp, R? D ) on we 
have that 

We may assume without loss of generality that C > 100. If \x— y\ > ^R 1 ~ 5e ° 
then 

g\x -y\> \CM >2M> \((x) - ((y)\. 
Therefore Lemma A. 2 implies, by possibly increasing C > 0, that 



CM 

(6.15) max|V5(6)|< 



R 

logR + 



Rl-eo 



CM 

< 



R l-6e D 



b£E* ' W ' _ R 1 ' 5 ^ 

on T~L. Trivially, 

M < max \h?(x) \ +max\h*(x)\ 

and by Lemma 3.3 we know that (EM p ) 1 / p = Oj^ p (R e ). This clearly gives 
(6.12). Combining (6.13) with (6.14) gives (6.11). ' □ 



44 



JASON MILLER 



6.3.3. Putting Everything Together. 

Proof of Theorem 1.2. To prove the theorem we just have to estimate 
the error terms from Proposition 6.2 on H. First of all, by Lemmas 3.3, 6.4 
we observe 

(6.16) (E[(£ yi H ]f /p = 0^ p (R 1+2 ^~ a ^) 
for a% = 1/10. Consequently, 

E[£ E l n ] =O x (i? 2 • R ic ^D-3 + R l+2e D -a iPB +c D e D -l^ 

(6.17) =0^{R CieD ~ aipB ) 

for ci < 100. Using Lemmas 3.3, 6.4 again, we see that 

E[£ e B l H ] =O x ( J R 1+2eD+7f + 2 - 1 )O x ( J R 1 / 2+CDeD - 1 + i^ +(7*-i/2)p B ) 

(6.18) =0 -( R cie D +2 le+2 -a lPB ^ 

the last line coming as /?b < 1- Finally, Lemma 6.4 clearly implies 

M 

E[£fl«] = 0_( i? l+72-7 fcjR e D + 7fc (l-p CD ) i? c 2 e D -l ) 
M 

(6.19) = ^ O x (ir 3£D ~ WCD ) 

fc=£+i 

for c 3 < 1000. The exponents from (6.17), (6.18), (6.19) are 

cie D - dips, cie D + 2j e+2 -aipB, c 3 e D - jipcD- 

Choosing I > 10 5 we see that the second and third exponents are negative 
and the first is clearly negative. □ 

We finish this section with the short proof of Theorem 1.3. 

Proof of Theorem 1.3. Assume that we still have the setup of the 
previous theorem except '0 = 0. Then there exists 5 > so that we can find 

a coupling of P^, , P^, so that h is harmonic in E on an event % (this is 
different from T~L in the proof of Theorem 1.2) with probability 1 — Oj^(R~ s ). 
Let g, h be the harmonic extensions of h, ~Eh^(x) from dE to E, respectively. 
For x £ E we have that 

Eh^(x) = Bh(x) = Bg(x)(l - l^ c ) + E7t(x)l,£ c 
= h(x) + E(h(x) - g(x))lfi c 
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since Eh^(x) = 0. Since ip S B^-(D), Lemma 3.3 implies both Eh 2 (x) = 

0^(R S / 2 ) and Eg 2 (x) = 0^(R 5 ^ 2 ). Consequently, an application of Cauchy- 
Schwarz yields 

E[\h(x)\ + \g(x)\]l nc =0^(R~ s /*), 
from which the theorem follows. □ 

7. The Central Limit Theorem. We will prove Theorem 1.1 in this 
section. The primary inputs are Theorem 1.2 and the main result of either 
[14] or [21]. Throughout, we let D C R 2 be a connected, bounded, smooth 
domain and for each n let D n = D n ^Z 2 . We will be dealing with both 
discrete and continuum derivatives, so to keep the notation consistent with 
the rest of the article we will still let V, A denote the discrete gradient and 
Laplacian, respectively, and use V and A for their continuum counterpart. 
We begin with a simple analysis lemma. Let /3 = (3(u) be as in the statement 
of Theorem 1.1. 

Lemma 7.1. For each e > there exists c > so that for all g £ C°°(D) 
we have that 

Y |V<7(6)| 2 < csup ||V 5 (x)|| 2 , Y IW)! 2 < c|| 
b&D* n xeD beD* 

Y l v s( & )l < c\\g\\ H 2+, {D) , Y \^9{x)\ < c\\g\\ H 3+e( D y 

bedD* xeD n 

Proof. The first claim is obvious since |V<7(6)| < - sup^g^ ||Vg(x)||. As 
for the second claim, we note that the Sobolev embedding theorem implies 
that for each e > there exists c, d > so that 

sup ||Vs(x)|| < c'\\Vg\\ H i+c {D) < c\\g\\ H 2+e, D \. 
xeD 

See, for example, from Proposition 1.3 in Chapter 4 of [27]. The final two 
claims are proved similarly. □ 

Let tp n (x) = nu-x. Fix a continuous function / : R 2 — > R and let h n have 
the law of the GL model on D n with h n (x) = f(x) + (p n (x) for x G dD n , 
v n,D _ y^n and for g e H 3+£( D } define 

& D (9)= £ a u (b)Vg(b)( V n ' D (b)-Vcp n (b)) 
beD* 

where a u (b) = E[V" \r]{b))} for ry ~ SEGGS„. Note that a u (b) depends only on 
the orientation of b. 
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Lemma 7.2. For each e > there exists c > so that for allg 6 H 3+e (D) 
we have 

C I 



Eexp(^' D (<7))<exp 



2 I II ||2 
oo + c \\9\\b 3 +z(D) 



Proof. Using summation by parts, we have 
|E# D (g)|< \E(h n (x)-Mx))^g(x)\+ci Yl \f( x b)Vg(b)\ 

x&D n b&dD* 



<C2\\j\\oo\\g\\H3+t(D)- 

In the final inequality we used that \~E(h n (x) — ip n (x))\ < ||/||oo) which is 
a consequence of Lemma 3.1, in addition to Lemma 7.1. The exponential 
Brascamp-Lieb inequality (Lemma 3.2) combined with (2.1) and the previ- 
ous lemma implies there exists C3, C4 > so that 



Bexp(^ D (g)) < exp 



< exp 



C3 



\9\\hs+*(d) + Yl ( V ^( 6 )) i 

6GD* 



2 1 i 1 1 2 

OO + \\9\\H3+e(D) 



Lemma 7.3. For each k > 4, the law of induces a tight sequence on 
H~ K (D) equipped with the weak topology. 

Proof. Fix k > 4. It suffices to show that for each 5 > there exists 
M = M(5) such that 

P[\K D \\h-^d)>M]<8 

as the Banach-Alaoglu theorem yields that the ball {Hgll/f^D) — M} is 
compact in the weak topology of H~ K (D). 

Let e = ^(k — 4) and k' = k — 1 — e>3. Let (<?&) be the eigenvectors of 
A on D, normalized to be orthonormal in L 2 (D), with negative eigenvalues 
(A/t) ordered to be non-increasing in k. Let = (1 — A) _K//2 ^. By (2.3), (gk) 



□ 



is an orthonormal basis of H K (D). As gk/(l — \k) <yl+e ^ 2 
we have that 

IklU'P) = 11(1" A)- (1+e)/2 5fc||^(D) 



;i-a)- 



1 



(l_ Afe )(l+e)/2- 

The Weyl formula implies that k/ (—Xk) tends to a constant c = cd depend- 
ing only on Z) as /c — > 00. Therefore there exists > cd so that 

\\k {l+e),2 9k\\ HK > {D) < c' D for all kN. 
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Combining the above with Chebychev's inequality yields 

n& D (9k)\ > M/k^W] = P[K' D (k (1+e)/2 9k)\ > Mk^] < exp(c-M^/ 4 ) 

where c = c(e, D, /). Consequently, letting A n M = n k {\Q D (g k )\ < M/fc( 1+e / 2 )/ 2 }, 
a union bound yields P^AjJf] — > 1 as M — > oo. Note that if g e H K (D), 
9 = Y, k a k9k for (a k ) £ I 2 , we have 

k 

since ff K (D)-convergence implies uniform convergence as k > 4 and is 
obviously continuous in the uniform topology on functions D n — > R. Let 

k 

On Aj^, note that 

1/2 / \ 1/2 



(D)MVN e . 



Consequently, 



sup \#%)\ >M^N e 

\\9\\h k (d)<^ 



This proves for every 5 > there exists M = M{5) sufficiently large so that 

P[\K' D Wh-Hd)>M]< 8 
for every n S N. □ 

Let r\ ~ SEGGS^, but thought of as a random gradient field on (^-Z 2 )*. Fix 
a base point x* £ dD and let x n be a point in dD n with minimal distance 
to x*. Set h n '°(x n ) = and let /i n '° be the function satisfying V/i n '° = rj. 
Let 

ftfo) =5> u (6)V<7(6)(7f (6) - Vp„(6)) for 5 G C^R 2 ). 

Corollary 2.2 of [ ! ] implies that for any gi,-..,g k G C^°(R 2 ) fixed, the 
random vector (£v(<?i)j • • • >£v(#*)) converges in distribution to a zero- mean 
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Gaussian vector (Zi,...,Zk) with covariance Cov(Zi,Zj) = (g%,gj)y for 
A = A(u, V) depending only on the tilt u and V. Note that our definition 
of differs from £ e in [14] in that we do not have a normalization of n _1 . 
The reason is that £ v operates on discrete gradients of C^°(R 2 ) functions, 
which themselves are of order n _1 . 

Lemma 7.4. There exists A depending only on V such that 

P[h n % Dn GB^(D n )] = l-0(n- s ). 

Proof. Let x e D n and Combining the exponential Brascamp- 

Lieb inequality with (4.6) yields for x G dD n that 

Eexp(/i n ' (x)-(^„(x)) = Eexp(h n '°(x)-h n '°(x n )-<p n (x)) < exp(Clogn) = n c . 
Assume without loss of generality that C > 1. By Chebychev's inequality, 

P[\h n '°(x) - tp n {x)\ > lOClogn] < n- g . 
Using a union bound we thus have 

P[ max \h n '°(x) - y n (x)\ > lOClogn] < 0(n" 8 ). 

x£dD n 

Consequently, taking A = IOC we have that P[h n >°\ dDn G B^(D n )\ = 1 - 
0(n- 8 ). □ 

By the same proof as Lemma 7.2 we have 
(7-1) Eexp(£ v (<?)) < exp(c\\g\\ 2 H3+e(D) ). 

PROOF of Theorem 1.1. Fix k > 4 and /: R 2 — > R continuous. Let 
h°, h,f be GFFs on D with respect to (•, -)y , A = A(u, V) as before, where 
h° has zero boundary conditions and the boundary condition of h* is given 
by f\dD- Let £y be a we&k-H~ K (D) subsequential limit of (£y D ). We will 
prove for any g±, ■ ■ ■ ,gt £ C°°{D) that 

(££(<?i), • • ■ ,£vGfc)) = ((/»', <&)4, ■ • • , (h f ,g k )$) 

since the continuity of £y implies that its law is determined by its projections 
onto C°°(D), a dense subset of H K {D). We will identify A at the end of the 
proof. To establish this, it suffices to show that 

(® D (9i), $ D {9k)) -4 ((/^ 51)4, • • • , (fc', <fc)4). 
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We will first prove the result for C^°(D), then using an approximation ar- 
gument generalize to C°°(D). 

By Lemma 7.4, with probability 1 — 0(n -8 ) we can apply Theorem 1.2 to 

,Pjq where ip = f + (p n and ^ = /i n '°|aD n . This implies the existence 
of e, 5 > independent of n such that we can couple h n , h n '° so that with 
h n the A^-harmonic extension of h n — h n '° from dD n (n~ e ) to D n (n~ e ), we 
have 

P[^ c ] = 0{n~ 5 ) where ?i = {h™ = h n in D„(™~ e )} 

for all n large enough. The reason that we see n~ e rather than n 1_e as in 
the statement of Theorem 1.2 is that D n = D n ^Z 2 , so all of our distances 
need to be scaled by a factor of rT x . Fix gi,...,gje G C^°(D) and assume 
that n is sufficiently large so that supp(gi), . . . , supp(gfc) C D n (n~ e ). On 
for each 1 < i < k we have that 

(7.2) £v' D (^) = &Gft) + 2 a«(6)V 5i (6)V^(6) = £ v (<7i)- 

fee A* 

The second equality follows from summation by parts and the A^-harmonicity 
of h n ; there is no boundary term since gi vanishes near dD n (n~ € ). 

Combining Lemma 7.2, (7.1), and the Cauchy-Schwarz inequality yields 

E|£v' D (^) = E|^ D ( 9i ) 

<[0(l)0(n- s )] 1/2 < 0{n- & / 2 ). 

Therefore (£y (51), . . • , £y (gk)) is a Gaussian vector with Cov(£y (#»), £v (#?')) = 
Cov((/t, 5i)y, (h,gj)Q) where fa is an j4-GFF on R 2 . Proposition 2.1 implies 
that /i restricted to D admits the decomposition h = h° + h where h° is 
a zero-boundary ^4-GFF on D and /i is a A -harmonic function. Integra- 
tion by parts implies that (h,gi)y = for all i, consequently the covariance 
structure of (£y(#i),-- .,£§(gk)) is the same as ((/i ,5i)^, . . . , (h°,g k )§). 

We now turn to the general case that g\, . . . , gk G C°°(D) do not neces- 
sarily have compact support in .D. Note that we can write 

gi = (gi -gi-di) + gi + g 

where gi is the A^-harmonic extension of gi from dD to D and gi G 
satisfies ||<7j — gi — gi(x)\\iji^ — Note that such an approximation exists 
since g%—gi G Hq(D). Since gi is harmonic with smooth boundary conditions, 
we have A@gi = o(l)n -2 uniformly in D. Thus summing by parts twice and 
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using that h n (x) = f(x)+tp n (x) on dD n , with f n the A^-harmonic extension 
of / from dD n to D n we have 

E£ v ' D (?i) = Y, a u(b)f(x b )VgiQ>) + = E au(b)Vf n (b)VMb) + o(l). 

Thus it is not difficult to see that if F denotes the A^-harmonic extension 
of / from 3D to D then 

lim E^(?i) = [ VFA u Vgi = [ VFA u V 9i 
n ^°° Jd Jd 

where A u is the diagonal matrix with entries j3 = (/Si,^). 
Applying summation by parts, 

bedD* n x&D n 

Note that the first summation on the right hand side is deterministic. Con- 
sequently, combining the Brascamp-Lieb inequalities with (2.1) implies 

Var(^' D (? 4 )) =0(1) E |A^(x)A^(y)|G n (x,y) 
= o(l) E ^ 2 E ^ 2 G„(x, 2/ ) = (l) 

where G n is the discrete Green's function on D n . This takes care of g{. We 
already know that the limiting behavior of (Jji), which leaves us to deal 
with gi — Tji —gi. Invoking Lemma 7.1 and the Brascamp-Lieb inequality, we 
have 



limsupE(£™' D (c/i -gi - gi)) 2 < Climsup E I V (#i _ 9i ~ 9i)Q>)\' 

beD* 

(7.3) =C\\gi-gi-gi\\ 2 Hl{D) <c5 2 . 



In the equality, we are using that gi — gi — gi £ C°°(D). We also have 

E|(/i°,3i -m -5i)vl 2 < c hi -9% -9i\\m(D) ^ cS i- 

Assume that (t$ D (gx), ■ ■ ■ , t$ D (gk)) and ((h°, 9i)y, ■ ■ ■ > (h°, g k )^) have 
been embedded into a common probability space so that lim n £y D Qi) = 
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(h°,gi)y almost surely for each 1 < i < k. By (7.3), 



Elimsup|£v'%i)-(/A3i)v- / VFA u V 9i 




=Elimsup |£ v ' (gt gi) - (h°,gi -gi~ gi)\ 



<25 



Since 5± > was arbitrary we therefore have that (£v(<?i)j • • • , £v(fl7c)) has 
the same distribution as 



We will now explain why A = A u , which will complete the proof. We will 
not spell out all of the details exactly in order to avoid repetition. Suppose 
that U C D is a smooth open subset, U n = U n ^Z 2 , and let ^ n be the 
function which is A@ harmonic in U n and is equal to h n (x) — (p n {x) in 
D n \ U n . Consider the auxiliary functional 



Exactly the same argument implies that £y' converges to a zero-boundary 
^4-GFF on U, say £y. Since ip n = £y D — £y ,C/ , as linear functionals, we also 
know that ip n has a limit, say tp = ^y — £y . It is not difficult to see that 

ip is A^ harmonic and depends on £y only through its values on D \U. 
Since £y is an ^4-GFF on D, we know that it admits the decomposition 

£y = £y + 1/) where ip is A -harmonic and £y is a zero-boundary yl-GFF on 
U independent of £y . Therefore we have that 



This implies that t/j = ip almost surely since a zero-boundary ^4-GFF plus 
an independent function does not have the law of a zero-boundary ^4-GFF. 




& U (9) = E *u(b)V(h n -<p n - iP n )(b)Vg(b). 



b£D* 



This finishes the proof of the theorem. 



□ 



APPENDIX A: DISCRETE HARMONIC FUNCTIONS 



Suppose that D C Z 2 is bounded and connected. We say that D satisfies 
the r-exterior ball property if for each x £ dD there exists y £ Z 2 such that 
x £ dB(y, r) and B(y, r) n D = 0. 
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Fig 4. The setup for the proof of Lemma A.l. Here, we choose a Mobius transforma- 
tion tp: C — > C determined by ip{z-\) — — 1, <p(zo) = 0, and tp(zi) — 1. In particular, 
ip(B(zo,s)) is mapped to D, the unit disk inC, and ip(B(za, s)\B) is sent to U, the upper 
half of the unit dtsk. 



Lemma A.l. Suppose that D C Z 2 is bounded, connected, and satisfies 
the r-exterior ball property. Let /?i,/32 > with Pi + P2 = 1/2. Let X be a 
random walk on D with Xq = x that jumps up and down with probability 
Pi and left and right with probability fi 2 - Let td = inf{t > : X t ^ D}, 
t s = inf{t > : \Xt — x\ = s}, and d = dist(x, dD). There exists a constant 
C = C(/3i,/3 2 ) > such that 

s A r 

Proof. We are first going to prove a related result for Brownian mo- 
tion, then explain how to deduce the corresponding result for random walk. 
Clearly, we may assume that d < s/4 and s < r/4. Suppose that z £ C, B is 
a ball of radius r, dist(z, B) = d, and zq is the point in dB closest to z. Let 
W be a Brownian motion initialized from w £ C, t^ = ini{t > : Wt 6 B}, 
and = inf{t > : | W t - z \ > s}. Finally, let u(w) = P w [t^ < rj]. We 
claim there exists C > independent of the setup such that if w E B(zq, s/4) 
then 

C 

u(w) < — dist(w,B). 

s 

We know that u solves the Dirichlet problem 

Au = in B(z ,s)\B, u|( 5B(2 . 0jS) )\ B = 1, n| (ajB )\ Bc(;ZOiS ) = 0. 

Let Z-i,Zi be the two points in dB n 8B(zq,s) and let 99 be the Mobius 
transformation satisfying ip(z-i) = — l,(p(zo) = 0,<f(zi) = 1. Let D = {z £ 
C : \z\ < 1}, U = {z e D : Im(z) > 0}, U x = {z 6 d\J : Im(z) > 0}, and 
U 2 = {z : lm(z) = 0}. Then 



<p(B(z , s) \B)=V, <p((dB(zo, s)) \B) = U U ip{{8B \ B c (z , s)) = U 2 . 
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Consequently, v = u o p~ solves 

At; = in U, = 1, v\u 2 = 0. 

Let Us = {z € dD : lm(z) < 0} and let v solve 

Av = in D, v\u 1 = 1, v\u 3 = — 1. 

By symmetry, v\u 2 = 0, hence v = v in U. Therefore u is the restriction of 
u = v o ip. Since 5 is harmonic in B(zo, s) with ||S|aB( ZOiS ) ||oo < 1 ft follows 
that u is Lipschitz in B(zq, s/2) with constant C/s > and C is independent 
of the setup. Fix w 6 -B(zo, s/4) and let wq be the point in dB closest to w. 
As w, wo £ B(zo, s/2), we have 

C C 
\u(w)\ = \u(w) — u(wo)\ < — \w — it;o| = — dist(ty, dB). 

This proves our claim, from which we will now deduce the lemma. Fix 
x £ D with dist(x, dD) = d. Since D satisfies the exterior ball property, 
there exists y £ 7? such that with B = B(y,r) we have B n D = and 
dist(i?,x) = d. By monotonicity, it suffices to show that P x [t s < tb] < 
Cd/(s A r), where tb = inf{£ > : X t G B}. This is slightly different than 
the setting for the Brownian motion case because r s is the first time X has 
distance s from x, whereas before we considered the first time W has distance 
s from zo, which was the point in B closest to z. By the obvious monotonicity 
of the problem, the desired bound in either case implies the corresponding 
bound in the other. The simple random walk estimate obviously follows 
from the Brownian motion estimate. The non-simple case follows since the 
Brownian motion estimate holds even after applying a non-degenerate linear 
transformation. □ 



Lemma A. 2. There exists a constant C > such that the following holds. 
Suppose that D C Z 2 is bounded, connected, and satisfies the r-exterior 
ball property. Suppose further that g: D — >■ R is a A' 3 -harmonic function, 
fiiifii > 0, such that there exists g > such that \g(x) — g(y)\ < kg if 
x, y G dD and \x — y\ = k. Then 



max | X7g (b) \ < Cg 



\ogR + 



R 



where R = diam(D). 
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Proof. Fix b = (x,y) G D*. Let Xt be a random walk, jumping up 
and down with rate j3i and left and right with rate 02, initialized from x, 
Y t = X t + {y - x), and a = inf{t > : X t G dD or Y t G dD}. Since 
g(Xt)-giYt) is a bounded martingale, the optional stopping theorem implies 
that 

\Vg(b)\ <E\g(X a )-g(Y a )\. 

Therefore it suffices to show that there exists C > so that for every b = 
(x, y) with x G dD and y G D we have 



|V<?(6)| < 



logi? + 



i? 



Let Yf be a simple random walk started from y, tjj = inf{t > : Yt G 3D}, 
and let p& = P^Y^i} — x\ > k]. By the previous lemma, we know that there 
exists C > so that < C /{k A r) when A; > 1. By summation by parts, 



R R-i 
I Vff(6)| < y](Pfc-i - <^Pk9< Cg 



k=i 



k=0 



R-1 



^ k 

k=l 



1 ^ 1 

^ ^ — ^ ^ 



k=r+l 



which proves the lemma. 

APPENDIX B: SYMMETRIC RANDOM WALKS 



□ 



Throughout we suppose that we have a continuous time random walk X t 
on Z 2 with time-dependent jump rates ct(b) satisfying 

< a < c t (b) < A < oo 

uniformly in the edges b and time t. Let p(s, t; x, y) be the transition kernel 
of X. One of the important tools in the analysis of such walks is the Nash- 
Aronson estimates, and their time dependent generalization, which give a 
comparison between p and the transition kernel of a standard random walk. 

Lemma B.l (Nash-Aronson Estimates). Let DCZ 2 and p(u,t; x,y) be 
the transition kernel of X, the random walk in D jumping with rates ct(b) 
stopped on its first exit from D. There exists C > 1,5 > depending only 
on a, A such that 



(B.l) 
(B.2) 



p(s,t;x,y) < 



p(s,t;x,y) > 



C 



1 V (t- 
5 



exp 



y\ 



C(l V (t-s)V2) 
for \x — y\ < y/t — s 



lV(t-s) 

provided s < t with \t — s\ < r 2 ~ 2e and x, y G B(xq, r) with B{xq, 2r) C D 
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Proof. In the time-independent setting with D = Z 2 these follow from 
the usual Nash-Aronson estimates, see [26]. The time-dependent extension, 
also for D = Z 2 , is proved in Propositions B3 and B4 of [14] and also 
Propositions 4.2 and 4.3 of [6]. 

This leaves us to handle the case that D ^ Z 2 . We can couple X, X 
together so that they are equal until X exits D. By the case for D = Z 2 , the 
probability that this happens before time s + r 2 ~ 2<E given X s = x £ B(xo,r) 
is of order exp(— r e /C) where C is the constant of that lemma. That is, 

sup sup \p(s, t; x, y) — p(s, t; x,y)\ = 0(exp(— r e /C)) for t < s + r 2 ~ 2<E , 

x£B(x ,r) yeZ 2 

which implies the result. □ 
Let Lt be the operator 



^ Lemma B.2 (Nash Continuity Estimate). Let D C Z 2 . Suppose that 
ft : [0, oo) x D — > R is a solution of the equation 



for all u < s < t with \t — u\ < r 2 2<E and x, y G B(xo, r) with B(xq, 2r) C D. 

Proof. In the same manner as Lemma B.l, this follows from the usual 
Nash continuity estimate for D = Z 2 , which is proved in the time indepen- 
dent setting in [26] and the time dependent version is Proposition B6 of 



The following is Proposition 4.1 of [6]. 

Lemma B.3 (Caccioppoli inequality). There exists C > depending only 
on a, A such that the following holds. If ft solves dtft = C-tft on [0, 2r 2 ] x 
B(xo,2r) then 



Ctf{x) = Y, c t(p)Vf(b). 




(B.3) \ft(x) - f.(y)\ 




[14]. 



□ 



(B.4) 
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We now prove a Beurling-type estimate for Xt. 

Lemma B.4. There exists constants ps, C > depending only a, A such 
that the following holds. Fix x G Z 2 ,r > and let H be a connected graph 
with Ht~)B(x, r) c ^ and dist(#, x) = d<r. Ifr r = inf{i > : \X t -x\ = r} 
and th = inf{i > : Xt G H} then 

P[r r <r H ]<C I- . 

The idea is to show that the probability that X runs around an annulus 
with inner and outer radii r, 2r, respectively, is strictly positive independent 
of r. The following auxiliary estimate will be useful for the proof. 

Lemma B.5. Fix a G (0, 1) and let F C B(x,aR). Let r = inf{i > : 
Xt ^ B(x, R)} and rp = inf{i > : Xt G F}. There exists C > depending 
only on a, A, a such that 

P x [tf < r] > C^J. 

Proof. Employing Lemma B.l in the final step, 

Px[tf < r] > P x [X t G F, r > t] = P x [X t G F] - P x [X t G F, r < t] 
> ^ 0(0, t; x, y) - E[E Xt [ P (t, t; X T ,y)]l {r < t} }) 



y&F 

>J2[ /»((>•/:•'•■//) E 



C / \X T -y\ 

ex P -^TTTTT^um 1 {^<*} 



lV(t-r) ^ V C(l V (t - r)V2) 



Obviously, |X T — y\ > (1 — a)i?. It is easy to see that x h-> Cx _1 exp(— (1 — 
a)i?C -1 z — 1 / 2 ) is maximized when x = (1 — a) 2 -R 2 /(4C 2 ). Consequently, if 
i = 7i? 2 for 7 > then 

_ r . v-^ / 5 C \ (5 C \ \F\ 

Px[TF ~ T] ~ 2-, " (1 - a) 2 i? 2 J " U"a^Fj^ 

where C" does not depend on a, 7 provided F C 5(a;, y/jR). This implies 
that there exists ao > such that the result holds whenever a G (0, ao). 
To see the general case, let x\, . . . ,x n be an ao/4-net of B(x, aR). If Xj G 
-6(0, ^p), then it follows from our a priori estimate that Xt hits B(xj, ^) 
before exiting B(x,R) with strictly positive probability po- Iterating this ar- 
gument, it follows that the probability that X t hits B(xk, ^p-) before exiting 
B(x,R) with strictly positive probability pi. The result is now follows by 
combining this fact once again with our a priori estimate. □ 
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Proof of Lemma B.5. The previous lemma implies the existence of 
pi > depending only on V such that the following holds. Let A(x, r\,r2) = 
{y : n < \x — y\ < V2} be the annulus with inner and outer radii r\,r<2 
satisfying T2 = 2ri. Let C be the event that X runs a full circle around 
A(x, T\,r2) after hitting dB(x, |rj) without hitting dA(x, r\, r%). Then 

(B.5) P[C] > pi > 0. 

Now set rk = 2 k d. The largest index m so that r m < r is [log2 a J • By (B.5), 
the probability that X makes it to distance r without running a full circle 
around one of the A(x, Vi, r^+i) is at most 

(1 - Pi) m < exp(l - \og 2 {r/d) Pl ) = C 

□ 
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